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ABSTRACT 

We have previously presented evidence for a varying fine-structure constant, a, in two inde- 
pendent samples of Keck/HIRES QSO absorption spectra. Here we present a detailed many- 
multiplet analysis of a third Keck/HIRES sample containing 78 absorption systems. We also 
re-analyse the previous samples, providing a total of 128 absorption systems over the red- 
shift range 0.2 < z a b s < 3.7. The results, with raw statistical errors, indicate a smaller 
weighted mean a in the absorption clouds: Aa/a = (—0.574 ± 0.102) x 10~ 5 . All three 
samples separately yield consistent and significant values of Aa/a. The analyses of low-z 
(i.e. z a bs < 1-8) and high-z systems rely on different ions and transitions with very different 
dependencies on a, yet they also give consistent results. We identify an additional source of 
random error in 22 high-z systems characterized by transitions with a large dynamic range 
in apparent optical depth. Increasing the statistical errors on Aa/a for these systems gives 
our fiducial result, a weighted mean Aa/a = (—0.543 ± 0.116) x 10~ 5 , representing 4.7 a 
evidence for a varying a. Assuming that Aa/a = at z a b s = 0, the data marginally pre- 
fer a linear increase in a with time rather than a constant offset from the laboratory value: 
a/a = (6.40 ± 1.35) x 10~ 16 yr _1 . The two-point correlation function for a is consistent 
with zero over 0.2-13 Gpc comoving scales and the angular distribution of Aa/a shows no 
significant dipolar anisotropy. We therefore have no evidence for spatial variations in Aa/a. 

We extend our previous searches for possible systematic errors, giving detailed analyses 
of potential kinematic effects, line blending, wavelength miscalibration, spectrograph tem- 
perature variations, atmospheric dispersion and isotopic/hyperfine structure effects. The latter 
two are potentially the most significant. However, overall, known systematic errors do not ex- 
plain the results. Future many-multiplet analyses of independent QSO spectra from different 
telescopes and spectrographs will provide a now crucial check on our Keck/HIRES results. 

Key words: atomic data - line: profiles - methods: laboratory - techniques: spectroscopic - 
quasars: absorption lines - ultraviolet: general 



1 INTRODUCTION 

High resolution spectroscopy of QSO absorption systems provides 
a precise probe of possible variations in several fundamental con- 
stants over cosmological distances and timescales. Dzuba, Flam- 
baum & Webb (1999a,b) and Webb et al. (1999) introduced and 
applied the new and highly sensitive many-multiplet (MM) method 
for constraining space-time variations of the fine-structure constant, 
a = e 2 /he. Using a MM analysis of 49 absorption systems from 
two independent samples of Keck/HIRES spectra, we recently re- 
ported statistical evidence for a variable a (Murphy et al. 2001a, 
hereafter MOla; Webb et al. 2001): Aa/a = (-0.72 ± 0.18) x 
10 , where Aa/a is defined as 



Aa/a = (a z - a )/a (1) 

for a z and ao the values of a in the absorption system(s) and in 
the laboratory respectively 1 . The absorption redshifts covered the 
range 0.5 < z a b s < 3.5. 

Here we significantly extend this work, presenting new mea- 
surements of Aa/a in a third independent Keck/HIRES sample 
of 78 absorption systems covering the range 0.2 < z a b s < 3.7. 
The present paper provides a detailed description of our methods, 
results and analysis of potential systematic effects. A companion 
paper, Webb et al. (in preparation), presents a summary of the main 
results. 

In the remainder of this section we outline the many multiplet 
technique of obtaining constraints on Aa/a from QSO spectra. We 
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1 a' 1 = 137.03599958(52) (Mohr & Taylor 2000). 
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describe the different QSO and atomic datasets in Section 2 and our 
analysis methods in Section 3. The new results, together with those 
from a re-analysis of the previous datasets, are presented in Section 
4. We consider both temporal and spatial variations in a. Section 
5 extends our search for systematic errors, concentrating on those 
difficult to quantify in Murphy et al. (2001b, hereafter MOlb). We 
summarize the present work in Section 6 and compare our new 
values of Aa/a with other constraints on a-variation. 



1.1 QSO absorption line constraints on Aa /a 

1.1.1 The alkali-doublet (AD) method 

The relative wavelength separation between the two transitions 
of an alkali doublet (AD) is proportional to a 2 (e.g. Bethe & 
Salpeter 1977). Savedoff (1956) first analysed AD separations seen 
in Seyfert galaxy emission spectra to obtain constraints on co- 
variation. Absorption lines from intervening clouds along the line 
of sight to QSOs are substantially narrower than intrinsic emission 
lines and therefore provide a more precise probe. Bahcall, Sargent 
& Schmidt (1967) first used AD absorption lines which seemed to 
be intrinsic to the QSO. Narrower still are absorption lines arising 
from gas clouds at significantly lower redshift than the background 
QSO. Wolfe, Brown & Roberts (1976) first analysed Mg II doublets 
from a damped Lyman-a system at z a bs ~ 0.524. Since then, sev- 
eral authors (e.g. Cowie & Songaila 1995; Varshalovich et al. 1996, 
2000) have applied the AD method to doublets of several different 
species (e.g. C IV, Si II, Si IV, Mg II and Al III). 

In Murphy et al. (2001c) we obtained the strongest current 
AD constraint on Aa/a by analysing 21 Si IV doublets observed 
towards 8 QSOs with Keck/HIRES. We found a weighted mean 



Aa/a = (-0.5 ±1.3) x 10" 



(2) 



The systematic error due to uncertainties in the Griesmann & Kling 
(2000) Si IV laboratory wavelengths is ±0.2 x 10~ 5 . Recently, 
Bahcall, Steinhardt & Schlegel (2003) analysed 44 [O III] dou- 
blet QSO emission lines from the Sloan Digital Sky Survey Early 
Data Release (Schneider et al. 2002). Their constraint, Aa/a — 
(—2 ± 1.2) x 10~ 4 , is an order of magnitude weaker than equation 
2. However, if one allows or expects spatial variations in a, such 
a simplistic comparison is difficult since the two results constrain 
Aa/a in very different environments (i.e. QSO emission regions 
and intervening absorption clouds). 



1.1.2 The many-multiplet (MM) method 

The many-multiplet (MM) method was introduced in Dzuba et al. 
(1999b) and Webb et al. (1999) as a generalisation of the AD 
method. The MM method constrains changes in a by utilizing 
many observed transitions from different multiplets and different 
ions associated with each QSO absorption system. The full details 
behind this technique are presented in Dzuba et al. (1999a, 2001b, 
2002). We presented a summary in M0 la and listed the many im- 
portant advantages the MM method holds over the AD method, in- 
cluding an effective order-of-magnitude precision gain. Therefore, 
we only outline the salient features of the method below. 

For small shifts in a (i.e. Aa/a <JC 1), the rest wavenumber, 
u> z , of a transition with a measured redshift, z, in a QSO spectrum, 
can be written as 
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Figure 1. The distribution of q coefficients in (rest) wavelength space for 
the low- and high-z samples. Note the simple arrangement for the low-z 
Mg/Fe II systems: the Mg transitions are used as anchors against which the 
large, positive shifts in the Fe II transitions can be measured. Compare this 
with the complex arrangement for the high-z systems: low-order distortions 
of the wavelength scale will have a varied and complex effect on Aa/a 
depending on which transitions are fitted in a given absorption system. In 
general, the complexity at high-z will yield more robust values of Aa/a. 
For the discussion in Sections 5.7.2 and 5.7.3 we define several different 'q- 
types' as follows: positive-shifters (q ^ 700 cm ), negative-shifters (q ^ 
-700 cm- 1 ), anchors < 300 cm" 1 ) and mediocre-shifters (300 ^ 
\q\ < 700 cm- 1 ). 



where ojq is the wavenumber measured in the laboratory 2 . If a z / 
Qo then 

-)'' I and y (—V 1 (4) 



aoJ 



are non-zero and the magnitude and sign of qi and qi determine 
the shift in the transition wavenumber. Since we only consider 
Aa/a -C 1, we may write 



LO z = LO + qxz 



(5) 



LO z = lo ± qiXz + qiVz 



(3) 



where we consolidate q\ and c/2 into q = q\ ±2c/2. The q coefficient 
represents all the relativistic corrections for the transition of interest 
and varies both in magnitude and sign from transition to transition. 
That is, if Aa/a / 0, the QSO absorption lines will be shifted in 
a distinct pattern with respect to their laboratory values. 

Fig. 1 illustrates the distribution of q coefficients for the dif- 
ferent transitions used in the MM method. The q coefficients are 
calculated in Dzuba et al. (1999a, 2001a,b, 2002) using many-body 
techniques to include all dominant relativistic effects. The values 
used in the present study are listed in Table 2. The uncertainties in 
q are typically < 30 cm -1 for the Mg, Si, Al and Zn transitions 
but are < 300 cm - for those of Cr, Fe and Ni due to the more 
complicated electronic configurations involved. However, it is im- 
portant to note that, in the absence of systematic effects in the QSO 
spectra, the form of equation 5 ensures that errors in the q coeffi- 
cients can not lead to a non-zero Aa/a. 



2 This equation does not strictly hold where level pseudocrossing occurs 
(Dzuba et al. 2001b, 2002). However, even in these cases, only the deriva- 
tive of uj with respect to a is required for Aa/a < 10~ 3 (i.e. for w z 



sufficiently close to loo). 
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Webb et al. (1999) first applied the MM method to 30 
Keck/HIRES QSO absorption systems (towards 14 QSOs) in the 
absorption redshift range 0.5 < z a b s < 1.6, fitting only the transi- 
tions of Mg I, Mg II and Fe II. The large positive q coefficients for 
the Fe II transitions are in sharp contrast to the small values for the 
Mg transitions (Fig. 1): the Mg lines act as anchors against which 
the large shifts in the Fell lines can be used to constrain Aa/a. 
The difference in q values (Aq fa 1200 cm -1 ) is much larger 
than for the Si IV transitions in the AD method (Aq fa 500 cm -1 , 
Murphy et al. 2001c) and so the MM method allows significantly 
increased precision. Additionally, analysing the 5 Fe II and 3 Mg 
transitions allows a statistical gain over the AD method. Webb 
et al. (1999) explicitly demonstrated this increased precision, ob- 
taining tentative evidence for a smaller a in the absorption clouds: 
Aa/a = (-1.09 ± 0.36) x 10~ 5 . 

In M0 la we analysed 18 independent Keck/HIRES absorp- 
tion systems in the range 1.8 < z a b s < 3.5. At these higher red- 
shifts, shorter (rest) wavelength transitions are shifted into the opti- 
cal bandpass. Fig. 1 shows the range of q coefficients available and 
their complex arrangement compared to the low-z Mg/Fell sys- 
tems. Such an arrangement should provide more robust estimates 
of Aa/a since any systematic effects in the QSO spectra will have 
a different effect on systems containing different transitions. De- 
spite this, the high-z: sample also yielded evidence for a smaller 
a in the past, consistent with the low-z: systems. Combining this 
with a re-analysis of the Webb et al. (1999) data and including 3 
low-z Mg/Fell systems from the new data, we found 4.1 a evi- 
dence for a variable a with 49 absorption systems over the range 
0.5 < z abs < 3.5: 

Aa/a = (-0.72 ± 0.18) x 10~ 5 . (6) 

Such a potentially fundamental result requires extreme 
scrutiny and warrants a thorough examination of possible system- 
atic errors. We carried out a detailed search for both instrumental 
and astrophysical systematic effects in MOlb but were unable to 
identify any which explained the result in equation 6. Upper limits 
were placed on potential effects from atmospheric dispersion and 
possible isotopic abundance evolution but we indicated that further 
exploration of these possibilities was required. 

In the present work we analyse a new, large sample of 
Keck/HIRES absorption systems using the MM method. We also 
re-analyse the previous datasets and significantly improve the anal- 
ysis of potential systematic errors. 



2 NEW AND PREVIOUS DATA 
2.1 Keck/HIRES QSO spectra 

The QSO observations were carried out by three independent 
groups, all using the high resolution spectrograph (HIRES, Vogt 
et al. 1994) on the Keck I 10-m telescope on Mauna Kea, Hawaii. 
We previously studied the first two samples in M0 la and we outline 
changes and additions to these samples in Section 2. 1 . 1 below. The 
new, large sample of 78 absorption systems is described in Section 
2.1.2. Combining all three samples gives a total of 128 absorption 
systems observed towards 68 different QSOs. 8 absorption systems 
(towards 8 QSOs) were observed by two different groups, bringing 
the total number of independent Keck/HIRES spectra to 76. The 
absorption redshifts cover the range 0.2 < z a b s < 3.7. 



2.1.1 Previous low- and high-z samples 

In MOla we analysed two independent Keck/HIRES datasets, 
termed the low- and high-redshift samples. The previous low-z 
sample was provided by C. W. Churchill and comprises 27 Mg/Fe II 
absorption systems in the spectra of 16 QSOs covering a redshift 
range 0.5 < z a b s < 1-8 (e.g. Churchill et al. 2000a). An additional 
Zabs = 1.4342 system was contributed by Lu et al. (1996) but these 
data have been re-reduced and are now included in the new sam- 
ple. The previous high-z: sample was provided by J. X. Prochaska 
& A. M. Wolfe and comprises 11 QSO spectra which contain 17 
damped Lyman-a systems (DLAs) in the range 1.8 < z a b s < 3.5 
and 3 lower redshift Mg/Fe II systems (Prochaska & Wolfe 1999). 
An additional DLA at z: abs = 2.625 was contributed by Outram, 
Chaffee & Carswell (1999). In the present work we add 2 further 
DLAs, again provided by J. X. Prochaska & A. M. Wolfe, bringing 
the total number of systems in the previous high-z sample to 23. 
Details of the data reduction for these previous samples are pro- 
vided in MOla. We have reanalysed these samples (Section 3) in 
a way consistent with the new, larger dataset described in the next 
section. 



2.1.2 The new sample 

The main improvement on our previous results is due to a new 
dataset kindly provided by W. L. W. Sargent. The observations were 
carried out during numerous observing runs from 1993 November 
to 1999 November. The spectral resolution for most spectra was 
FWHM fa 6.6 km s" 1 (R = 45000). All spectra were reduced 
with the MAKEE 3 package written by T. Barlow. QSO exposure 
times ranged from 1000-6000 s depending on the QSO magnitude 
and between 3-16 separate exposures were combined to form the 
final spectrum. 

Thorium-argon (ThAr) emission lamp spectra, taken both be- 
fore and after the QSO exposures, were used to calibrate the wave- 
length scale. MAKEE uses an internal list of ThAr air wavelengths 
derived using the Edlen (1966) formula from the measured vacuum 
thorium atlas of Palmer & Engleman (1983) and the argon atlas of 
Norlen (1973). MAKEE then converts the final wavelength scale to 
vacuum wavelengths using the Cauchy formula from Weast (1979) 
and not the (inverse of the) Edlen formula. As we noted in MOlb, 
this distortion of the wavelength scale could lead to systematically 
non-zero values of Aa/a. For example, the effect on a z a b s = 1.0 
Mg/Fe II system would be a spurious Aa/a fa —0.25 x 10~ 5 . 
We have corrected the wavelength scales of the reduced spectra to 
avoid such a distortion (see MOlb for details). We test the accuracy 
of this wavelength scale in Section 5.3. The effect of slight miscal- 
ibrations (if any are present) on the overall Aa/a values is found 
to be negligible. 

For each QSO spectrum we identified all absorption systems 
which contained enough MM transitions to yield meaningful con- 
straints on Aa/a. We selected ~500-kms _1 sections around each 
transition and defined a continuum by fitting a Legendre polyno- 
mial. 1 a error arrays were generated assuming Poisson counting 
statistics. The average S/N per pixel for these spectral regions of 
interest ranged from 4-240, with most sections having S/N ~ 30 
(see Fig. 2). 

The new sample comprises 78 absorption systems, observed 
towards 46 QSOs, covering a wide redshift range, 0.2 < z a b s < 

3 See http://www2.keck.hawaii.edu:3636/inst/hires/makeewww. 
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Figure 2. The distribution of S/N and central wavelength for the spectral 
regions of interest in the previous low-z, previous high-z and new samples. 

3.7. Like the previous samples, the new sample conveniently di- 
vides into low-z (z abs < 1.8) and high-z (z abs > 1.8) subsamples. 
Table 1 shows the number of times each transition is used in the dif- 
ferent subsamples. Mg/Fe II systems dominate the low-z region: 35 
out of the 44 systems below z abs = 1.8 contain only Fe II and Mg I 
or II, and, out of the remaining 9, only 1 system does not contain 
transitions of these ions. We provide an example Mg/Fe II absorp- 
tion system, together with our Voigt profile fit, in Fig. 3. Note that 
fitting a large number of transitions allows the velocity structure 
of the absorption system to be reliably determined. This is further 
facilitated by the range in Fell oscillator strengths and the pres- 
ence of the weak Mg I line. The high-z systems are characterized 
by a diverse range of transitions from different ionic species. Table 
1 shows that the strong Sill A1526, Al II A1670 and Fell A1608 
transitions are the most common. Fig. 4 illustrates a particularly 
high column density absorber which contains the lower abundance 
species, Cr II, Ni II and Zn II. Theoretically, the large magnitude of 
the q coefficients for the transitions of these species makes them 
very important for constraining Aa/a (see Fig. 1). However, since 
their optical depths are typically < 0.3 and since they are not de- 
tected in all high-z systems, their influence over Aa/a is reduced. 

2.2 Atomic data 

We summarize all relevant atomic data for the MM transitions of 
interest in Table 2. This is an updated version of table 1 in M0 la, 
the main changes being the following: 

(i) q coefficients. As noted in Section 1.1, the qi and q 2 coeffi- 
cients used in previous works are now replaced by a single value, 
q = qi + 2q2. The values of q for Mg I and Mg II are taken from 
Dzuba et al. (1999a,b) and those for Al II and Al III are from Dzuba 
et al. (2001b). All other q coefficients are from updated calcula- 
tions of the type detailed in Dzuba et al. (2002). Conservative un- 
certainties are estimated by comparing other calculated quantities 



Table 1. The frequency of occurrence for each transition in the different 
samples of QSO absorbers. For the new sample, we define low-z as z abs < 
1.8 and high-z as z abs > 1.8. The previous low-z sample adheres to this 
definition but the previous high-z sample does contain 3 low-z systems. 
The low-z samples are dominated by the Mg/Fe II transitions. The high-z 
systems contain a diverse range of transitions and species but the strong Si II 
A1526, Al II A1670 and Fe II A1608 transitions are the most common. 
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(e.g. energy intervals, g-factors) with their experimental values and 
by comparing the results of several different calculation techniques. 
Thus, the errors quoted in Table 2 should be reliable but should not 
be treated as statistical. 

Of particular note are large changes to the q coefficients for 
the Fell A1608 and Nill lines. The q coefficient used in MOla 
for Fell A1608 was q = 1002 cm -1 but this was based on an 
incorrect excited state electronic configuration reported in Moore 
(1971), 3d 6 4p y 6 Py 2 - The configuration shown in Table 2 is taken 
from Dzuba et al. (2002) who also calculate the new q coefficient, 
q = —1200cm -1 . Since Fell A1608 is prominent in the high-z 
systems (see Table 1), this change has an important effect on the 
values of Aa/a in the previous high-z sample. We discuss this in 
detail in Section 4. The smaller changes to the Ni II q coefficients 
are due to the complicated electronic configuration of Ni II (Dzuba 
et al. 2001b, 2002). In this case, however, the effect on the values 
of Aa/a in the previous high-z sample is very small since the Ni II 
transitions are weak and only appear in ~25 per cent of our high-z 
systems (see Table 1). 

(ii) Oscillator strengths. In order that the parameters derived 
from our Voigt profile fits to the QSO data can be compared with 
those of other studies, we have used the set of oscillator strengths 
recommended in the DLA database 4 of Prochaska et al. (2001) (and 
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Figure 3. Mg/Fe II absorption system towards Q1437+3007 at z^bs = 1.2259. The data have been normalized by a fit to the continuum and plotted as a 
histogram. Our Voigt profile fit (solid curve) and the residuals (i.e. [data] — [fit]), normalized to the lcr errors (horizontal solid lines), are also shown. The 
tick-marks above the continuum indicate individual velocity components. Note the range of linestrengths in Fe II, facilitating determination of the velocity 
structure. The large number of Fe II transitions and the large number of velocity components allows for tight constraints to be placed on Ace/ce. The Mg I 
A2852 and Fe II A2374 transitions allow Aa /a to be constrained by the velocity m km s — 1 components whereas the Mg II and remaining Fe II transitions 
constrain An/a with the velocity fa —35 km s component. 



references therein). These changes have a negligible effect on our 
values of Act/ a. 

(iii) Aim hyperfine structure. In M0 la we used composite 
wavelengths to fit the Al III profiles. However, Griesmann & Kling 
(2000) clearly resolved the hyperfine components in their labo- 
ratory Fourier transform spectra (FTS) and found that the hyper- 
fine splitting is ~ 2.5 km s -1 , comparable with the resolution of 
Keck/HIRES. In MOlb we found that the composite values were 
adequate for our previous high-z sample, primarily because Al III 
only appeared in 6 absorption systems. However, we now include 
the Al III hyperfine components listed in Table 2 in our Voigt profile 
fits to the previous and new samples. 



Aside from Al, all the elements in Table 2 have several natu- 
rally occurring isotopes and so each QSO absorption line will be a 
blend of absorption lines from each isotopic component. However, 
the isotopic structures are only known for the transitions of Mg I 
(Hallstadius 1979) and Mgll (Drullinger, Wineland & Bergquist 
1980). The latter authors also determined the hyperfine splitting 
constant for the Mg II transitions. Pickering et al. (1998) used labo- 
ratory FTS to derive the wavelengths of the individual isotopes and 
main hyperfine components for these transitions and these are the 
values given in Table 2. 

In M0 la we estimated the isotopic structures for the two Si II 
transitions by scaling the isotopic shifts for Mg II A2796 by the 
mass shift, 

AtOi oc luq /m 2 , (7) 



4 Available at http://kingpin.ucsd.edu/~hiresdla. 



for Auii the shift in wavenumber for isotope i where m is the 
atomic mass. We give these estimated isotopic structures in Table 



6 M.T. Murphy, J. K. Webb &V.V. Flambaum 



Q0528-2505, z e =2.77, z. =2.1406 

' em ' abs 



p-i — i — r 



-i — i — i — r 



-i — i — i — r - 



~ I 1 1 1 1 1 1 1 1~ 



T~TF 



-i — i — i — r - 



t — i — i — r 



ID 

d 
o 



-I I II I I II Mill III It 







:Sill1526 



±Nill1741 

"i " i " " i " 



::Crll2062 



I I I I | I I I I | 4 



++++ 



I I I I 



| I I I I 



X 
Z5 



in 
d 

o 



CD 
N 



jji Jin n _n ,n 
' lni "imiiu' ir" v r : 



II II I INI I I I III 



I 



« « « !1 a n_n A _ _ _n „ II 1 1 -nrniJ 1 ! II _ 



I III 



:Fell1 608 



::Nill1751 



::Crll2066 



i i i 



i i i i 1 i i i i | 




i i | i i i i | 



O 

E ^ 

O LO 

Z d 
o 



LO 

d 
o 



III I II I Nil I 1 1 III " 



lr- B " , --f 



== - AJ I I I II I INI I I 



tAIII1670 W 



::Sill 1 808 



::Fell2344 ^ 



i i i 



■4—1- 



I I I I I I I I I 



i i i i " | i i i i 



n. ri. 



iiL n n. n.n 

-l^r-y-Lr"; 



II I I 



L w ^ n . .__n 1 

r it ir \J"i'\|LiHVJIi| 1 JVi lV -.'!!l W 



tNill1709 

. ■ • r • 




±Crll2056 



v 



f Fell2382 L / ■ 



100 100 -100 100 -100 

Velocity (km/s) [z = 2.141050] 



100 



Figure 4. Heavy element absorption lines in the damped Lyman-a system towards Q0528— 2505 at z a b s = 2.1406. The data have been normalized by a fit to 
the continuum and plotted as a histogram. Our Voigt profile fit (solid curve) and the residuals (i.e. [data] — [fit]), normalized to the 1<t errors (horizontal solid 
lines), are also shown. The tick-marks above the continuum indicate individual velocity components. Note the large range in linestrengths for the different 
transitions. The velocity structure is well-constrained by Sin A1526, Alii A1670 and Fell A1608 while the Nill and Cm transitions are strong enough to 
constrain Aa/o with the velocity km s components. This system is at low enough redshift for the Fe II AA2344 and 2382 transitions to be detected. 
These combine with Si II A 1526, Al II A 1670 and Fe II A 1608 to constrain Aa/a with the components at velocity 80 km s . 



2. Our approximation may break down if the specific isotopic shifts 
for the Si II transitions are very different to those of Mg II A2796. 
Very recent theoretical calculations (Berengut, Dzuba & Flambaum 
2003) suggest that these are of opposite sign to the normal mass 
shift and so the shifts in Table 2 are probably overestimates. We 
discuss the effect this may have on our values of Aa/a in Section 
5.5. 

At the time of our calculations, no measurements or theoret- 
ical estimates of the isotopic structures existed for the other tran- 
sitions (to our knowledge) and so we used the composite values 
listed in Table 2. Since that time, measurements (Matsubara et al. 
2002) and theoretical estimates (Berengut et al. 2003) of the Zn II 
isotopic structure have become available. In Section 5.5 we discuss 
the potential effects of differential isotopic saturation and strong 
cosmological evolution of the isotopic abundances on our values of 
Aa/a. 



Finally, we note that the Whaling et al. (1995) Ar II wavenum- 
bers used to calibrate the Sin, AllI and Aim wavelengths are 
systematically larger than those of Noiien (1973) used to cali- 
brate the wavelengths of the remaining lines. The difference be- 
tween the calibration scales is proportional to the wavenumber: 
Slo — 7 x 10~ 8 oj. In MOlb we note that calibration errors of 
this type will be absorbed by the redshift parameters of the QSO 
absorption lines and so, as long as all laboratory wavelengths are 
normalized to the one calibration scale, values of Aa/a will be un- 
affected. We have therefore normalized all laboratory wavelengths 
in Table 2 to the Norlen (1973) calibration scale. 

3 ANALYSIS 

We have detailed our analysis procedure in M0 la and so only pro- 
vide a short summary here. For each absorption system we con- 
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Table 2. Atomic data for the MM transitions in our analysis. Information for isotopic and hyperfine components is given in italics. Column 2 shows the 
mass number, A, for each species. The origin of the laboratory wavenumbers (ujq) and wavelengths (Aq) is summarized in MO la. Columns 5 and 6 show 
the updated ground and excited state electronic configurations. The ID letters in column 7 are used in Table 3 to indicate the transitions used in our fits to 
each absorption system. Column 8 shows the ionization potential for the relevant ion, IP + , and for the ion with a unit lower charge, IP - . Column 9 shows 
the oscillator strengths, /, from the DLA database 4 of Prochaska et al. (2001) or the relative strengths of the isotopic (Rosman & Taylor 1998) or hyperfine 
components (italics). The q coefficients are from Dzuba et al. (1999a,b, 2001b, 2002) and the uncertainties are discussed in the text. The Sill, Alii and 
Al III wavenumbers have been scaled from their literature values due to the Norlen/Whaling et al. calibration difference. 
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Pickering, Thorne & Webb (1998); "Griesmann & Kling (2000); dickering et al. (2000); "Pickering et al. (2002); e Nave et al. (1991). 
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struct multiple velocity component Voigt profile fits to the data us- 
ing VPFIT (v5) 5 . Each velocity component is described by three pa- 
rameters: the column density, the Doppler width or b parameter and 
the redshift, z abs , of the absorbing gas. We have modified VPFIT to 
include Aa/a as a free parameter: for each transition we alter the 
rest wavelength using equation 5 so that all velocity components 
shift in concert. To constrain Aa/a one must tie together the z a b s 
parameters of the corresponding velocity components in different 
transitions. This assumes negligible proper motion between the ab- 
sorbing gas of different ionic species. We discuss this assumption 
in detail in Section 5.1. 

We minimize \ 2 simultaneously for all fitted transitions to 
find the best fitting value of Aa/a. Parameter errors are calculated 
from the diagonal terms of the final parameter covariance matrix 
(Fisher 1958). We have performed Monte Carlo simulations of ab- 
sorption systems of varying complexity (i.e. different numbers of 
velocity components and transitions fitted) and with a range of in- 
put values of Aa/a to verify that VPFIT returns the correct values 
and 1 a errors (see also Section 5.7). 

We impose two self-consistency checks on each absorption 
system before accepting a value of Aa/a. Firstly, we require the 
value of x 2 per degree of freedom, to be ~1. Secondly, we 
construct two fits to the data with different constraints on the b 
parameters: (i) entirely thermal broadening and (ii) entirely tur- 
bulent broadening 6 . We require that the values of Aa/a derived 
from both fits be consistent with each other since the choice above 
should not greatly affect Aa/a (i.e. the measured position of ve- 
locity components should not change systematically). We rejected 
only 4 systems in this way. 



4 RESULTS 

4.1 Raw QSO results 

We present the \ 2 minimization results in Table 3 for the previous 
and new samples separately. We list the QSO (B1950) name, the 
emission redshift, z cm , the nominal absorption redshift, z a bs, and 
Aa/a for each absorption system with the associated 1 a statisti- 
cal error. The transitions fitted in each system are indicated with 
the ID letters defined in Table 2. We present basic statistics for the 
different samples in Table 4: (Aa/a) w is the weighted mean with 
1 a error, {Aa/a) is the unweighted mean, So is the significance 
of the departure of the weighted mean from zero, rms is the root- 
mean-square deviation from the mean of Aa/a, (S(Aa/a)) is the 
mean 1 a error and xt is the value of \ 2 P er degree of freedom 
with {Aa/a) w as the model. 

The statistics in Table 4 indicate a smaller a over the redshift 
range 0.2 < z a bs < 3.7 at the 5.6 a significance level: 

Aa/a = (-0.574 ±0.102) x 1(T 5 . (8) 

We note the consistency between the previous low-z, previous 
high-z and new samples. Breaking the new sample down into low- 
z (zabs < 1-8) and high-z (z a b s > 1-8) subsamples also yields 

5 Available at http://www.ast.cam.ac.uk/~rfc/vpfit.html 

6 This is a minor departure from our previous analysis. In MO la we con- 
structed a third fit for each absorption system where the thermal and tur- 
bulent component of the b parameters were determined by goodness of fit. 
Such a fit is more computationally intensive and, in general, the thermal and 
turbulent b parameters are poorly determined. In MO la we found that this 
fit did not lead to any inconsistent values of Aa/a and so fits (i) and (ii) 
provide a sufficient consistency check. 



Table 3. The raw results from the \ 2 minimization procedure. For each 
absorption system we list the QSO emission redshift, z em , the nominal 
absorption system redshift, z abs , the transitions fitted and the value of 
Aa/a with associated 1 a statistical error. A '*' after Aa/a indicates sys- 
tems included in the 'high-contrast' sample. An additional random error of 
2.09 X 10~ 5 should be added in quadrature to these systems to form the 
fiducial sample (see Section 4.2). 
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Table 3 - continued. The new sample. 
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Table 3 - continued. The new sample (continued). 
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2.551 + 5.392 


1337+1121 


2.97 


2.7955 


dhl 


4.103 + 8.538 


1425+6039 


3.20 


2.7698 


dgl 


-0.688 ± 1.843 






2.8268 


dgl 


0.433 + 0.827 


1437+3007 


3.00 


1.2259 


abcnopqr 


0.308 ± 1.460 


1442+2931 


2.76 


2.4389 


dgl 


-0.882 ± 1.473 


1549+1919 


2.83 


1.1425 


cfnpqr 


-0.076 + 0.671 






1.3422 


nopq 


-0.740 + 1.232 






1.8024 


def 


-3.050 + 2.473 


1626+6433 


2.32 


0.58596 


bcr 


-1.977 + 4.529 






2.1102 


defgl 


-0.705 ± 1.068 


1634+7037 


1.34 


0.99010 


abcdnopqr 


-2.194 + 1.343 


1850+4015 


2.12 


1.9900 


ghijklmnopqrvw 


-1.663 ± 0.859* 


1946+7658 


3.02 


1.7385 


efhijmntu 


-0.212 ± 1.857 






2.8433 


dgl 


-4.743 ± 1.289 



Object 


^em 


■^abs 


Transitions" 


A / /-in — 5 \ 

Aa/a (10 ) 


2145+0643 


1.00 


0.79026 


bcnopqr 


0.087 ± 0.589 


2206-1958 


2.56 


1.0172 


abcqr 


1.354 ± 0.883 






2.0762 


dl 


1.429 + 3.022 


2231-0015 


3.02 


1.2128 


abcopqr 


1.223 ± 1.465 






2.0653 


efghijklmstuvw 


1.707 + 1.249* 






2.6532 


dglstu 


-3.348 ± 1.904* 


2233+1310 


3.30 


2.5480 


dgl 


2.942 + 5.207 






2.5548 


defgl 


1.732 ± 6.349 






3.1513 


dgl 


-4.005 + 3.301 


2343+1232 


2.52 


0.73117 


abcoqr 


-1.211 + 0.975 






1.5899 


acdefgqr 


0.453 ± 1.187 






2.1711 


cdefnp 


-0.961 ± 1.295 






2.4300 


dgijklnqrstuvw 


-1.224 + 0.389* 


2344+1228 


2.77 


1.0465 


abcinopqr 


-0.747 + 1.530 






1.1161 


abcopqr 


0.009 ± 1.963 






2.5378 


dgl 


-3.205 + 2.094 



"Transitions identified as in Table 2; b This absorber contributed by Outram 
etal. (1999). 
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Figure 5. Histograms of raw results from Table 3. The distribution of 
Aa/a for the previous low-z (dot-dashed line), previous high-z (dotted 
line) and new (dashed line) samples are shown together with the overall 
distribution (solid line, binned differently in Fig. 15a). 



consistent results. We also see an overall agreement between the 
weighted and unweighted means, indicating that our results are not 
dominated by a small number of highly significant points. This is 
supported by the histogram in Fig. 5 which shows a symmetric, 
roughly Gaussian distribution for the values of Aa/a. 

We illustrate the distribution of Aa/a over redshift and cos- 
mological time in Fig. 6. The upper panel shows the raw values of 
Aa / a (with their 1 a error bars) for each sample separately. The 
middle panel shows an arbitrary binning of each sample such that 
each bin has an equal number of points (aside from the highest z 
bin which also contains the remaining points after dividing the to- 
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Table 4. Statistics for the raw QSO results and tests for systematic errors. The previous samples are described in Section 2.1.1. For the new and total 
samples, 'low-2' means 2 abs < 1.8 and 'high-2' means 2 abs > 1.8. For each sample, (z abs ) is the mean absorption redshift, iV abB is the number of 
absorption systems, (Aa/«) w is the weighted mean, (Aa/a) is the unweighted mean, So is the significance of the weighted mean (with respect to zero), 
rms is the root-mean-square deviation from the mean, {5 (Aa/a)) is the mean 1 a error and \l is X 2 per degree of freedom using (Aa/a) w as the model. 



Sample 


( z abs) 


^abs 


<Aa/a) w 


(Aa/a) 


So 


rms 


(5(Aa/a)) 


xl 


QSO results 


















Raw results (Fig. 6) 


















Prev. low-2 


1.00 


27 


-0.513 ±0.224 


-0.305 ± 0.304 


2.28(7 


1.579 


1.554 ±0.153 


0.981 


Prev. high- 2 


2.17 


23 


-0.672 ± 0.244 


-0.656 ± 0.601 


2.76(7 


2.880 


1.951 ±0.232 


3.082 


New low-z 


1.11 


44 


-0.537 ±0.159 


-0.652 ± 0.321 


3.37cr 


2.131 


1.757 ± 0.204 


1.123 


New high-2 


2.60 


34 


-0.623 ± 0.224 


-0.053 ± 0.473 


2.78(7 


2.755 


2.604 ±0.318 


1.533 


New total 


1.76 


78 


-0.566 ±0.130 


-0.391 ± 0.276 


4.36<7 


2.441 


2.126 ±0.186 


1.284 


low-2 total 


1.07 


74 


-0.539 ±0.130 


-0.562 ± 0.226 


4.24cr 


1.941 


1.663 ±0.134 


1.051 


high-2 total 


2.50 


54 


-0.637 ±0.172 


-0.227 ± 0.389 


3.72cr 


2.862 


2.400 ± 0.224 


2.179 


Raw sample 


1.67 


128 


-0.574 ±0.102 


-0.421 ±0.210 


5.62 a- 


2.379 


1.974 ±0.127 


1.515 


Removing extra scatter at high- 


z (Figs. 7 and 8) 














High-contrast sample removed 


2.65 


32 


-0.518 ± 0.347 


-0.011 ± 0.517 


1.49(7 


2.926 


3.070 ± 0.308 


1.207 


High-contrast errors increased 


2.50 


54 


-0.563 ± 0.291 


—0.227 ± 0.389 


1.94(7 


2.862 


2.881 ± 0.191 


1.103 


All high-2 errors increased 


2.50 


54 


-0.601 ±0.328 


-0.227 ±0.389 


1.83(7 


2.862 


2.986 ±0.194 


1.000 


High-2 (j-clip 


2.54 


44 


-0.489 ±0.191 


0.273 ±0.341 


2.57cr 


2.261 


2.559 ±0.263 


0.997 


Fiducial sample 


1.67 


128 


-0.543 ±0.116 


-0.421 ±0.210 


4.66er 


2.379 


2.177 ±0.124 


1.065 


Systematic error tests 


















ThAr calibration test: (Aa/a) 


ThAr (Fig. 


15) 














Prev. low-2 


1.00 


27 


0.016 ±0.021 


0.034 ± 0.045 


0.74cr 


0.234 


0.129 ±0.008 


3.697 


Prev. high-2 


2.10 


18 


-0.013 ± 0.020 


-0.015 ± 0.059 


0.66ct 


0.250 


0.165 ±0.043 


3.212 


New low-z 


1.14 


41 


-0.008 ±0.014 


-0.003 ± 0.036 


0.59 cr 


0.229 


0.123 ±0.016 


3.095 


New high-2 


2.63 


32 


U.U19 ± 0.U14 


0.U16 ± U.U52 


1.38 cr 


0.293 


n 1 jo I rt r\ o o 

U.14o ± [J.U22 


2.714 


New total 


1.79 


73 


0.006 ± 0.010 


0.005 ±0.030 


0.57 cr 


0.259 


0.134 ±0.013 


2.915 


low-2 total 


1.08 


71 


-0.006 ±0.011 


U.UUo ± U.Uz i 


U. J4 (7 




U. 1Z4 ± U.U1U 


3. ZOO 


high-2 total 


2.52 


47 


0.014 ±0.012 


v_j.uj.iJ _i_ yj.vj'-ti. 


1.24(7 


0.284 


\J.±tJZ? _1_ U.Uii 


2.827 


Total 


1.66 


118 


0.004 ± 0.008 


0.009 ±0.023 


0.48 cr 


0.253 


0.138 ±0.011 


3.091 


Composite wavelengths only: (A«/«) comp (Fig. 


16, top panel) 












Prev. low-2 


1.00 


27 


-0.802 ±0.227 


-0.557 ±0.334 


3.54(7 


1.736 


1.567 ±0.156 


1.094 


Prev. high-2 


2.17 


23 


-0.736 ±0.243 


-0.731 ± 0.608 


3.03(7 


2.918 


1.929 ± 0.230 


3.178 


New low-2 


1.11 


44 


-0.550 ±0.161 


-0.609 ± 0.329 


3.42cr 


2.184 


1.779 ±0.209 


1.060 


New high-2 


2.60 


34 


-0.650 ±0.224 


n nor 1 j~i i r/ 1 

—0.085 ± 0.456 


2.90cr 


2.660 


2.636 ± 0.326 


1.457 


New total 


1.76 


78 


-0.584 ±0.131 


-0.380 ±0.274 


4.47cr 


2.417 


2.152 ±0.191 


1.218 


low-2 total 


1.07 


74 


-0.650 ±0.128 


— U.04o ± Ki.Zoo 


J.U5 (7 


Z.U4J 


1 fiQ1 _l_ fl 1 QQ 

l.ool ± U.loo 


t.uyo 


high-2 total 


2.50 


54 


-0.667 ±0.172 


—0.251 ± 0.380 


3.89(7 


2.792 


2.411 ± 0.229 


2.129 


Total 


1.67 


128 


-0.656 ±0.103 


-0.481 ±0.212 


6.39 er 


2.397 


1.989 ±0.129 


1.520 


Strong isotopes only: (Aa/a) 


so (Fig. 16, 


middle panel) 












Prev. low-2 


1.00 


27 


-1.204 ±0.228 


-0.885 ± 0.353 


5.29(7 


1.834 


1.567 ±0.157 


1.257 


Prev. high-2 


2.17 


23 


-0.738 ±0.244 


-0.734 ± 0.619 


3.03(7 


2.967 


1.928 ± 0.226 


3.210 


New low-2 


1.11 


44 


-0.860 ±0.159 


-0.834 ±0.323 


5.40(7 


2.145 


1.768 ± 0.206 


1.279 


New high-2 


2.60 


34 


-0.634 ±0.224 


— U.lw / it U.4DO 


Z.64 cr 


L. /UU 


Z.OZO ± U.ozo 


1 AQA 


New total 


1.76 


78 


-0.784 ±0.130 


-0.504 ±0.275 


6.04(7 


2.432 


2.142 ±0.190 


1.363 


low-2 total 


1.07 


74 


-0.984 ±0.128 


-0.913 ± 0.239 


7.71cr 


2.055 


1.674 ±0.137 


1.300 


high-2 total 


2.50 


54 


-0.638 ±0.171 


-0.231 ± 0.384 


3.72cr 


2.824 


2.404 ±0.228 


2.143 


Total 


1.67 


128 


-0.861 ±0.102 


-0.626 ±0.215 


8.41 a 


2.433 


1.982 ±0.128 


1.663 


Atmospheric dispersion corrected: (Aa/a 


)adc (Kg. 17) 












Prev. low-2 


1.00 


27 


-0.098 ±0.225 


0.142 ±0.339 


0.43 cr 


1.762 


1.554 ±0.153 


1.208 


Prev. high-2 


2.17 


23 


-0.637 ±0.244 


-0.668 ± 0.583 


2.61(7 


2.794 


1.951 ±0.232 


2.979 


New low-z 


1.11 


44 


-0.416 ±0.159 


-0.557 ±0.319 


2.61cr 


2.118 


1.757 ±0.204 


1.050 


New high-z 


2.60 


34 


-0.740 ±0.224 


-0.189 ± 0.472 


3.31cr 


2.755 


2.604 ±0.318 


1.545 


New total 


1.76 


78 


-0.525 ±0.130 


-0.397 ± 0.274 


4.48cr 


2.423 


2.126 ±0.186 


1.266 


low-z total 


1.07 


74 


-0.288 ±0.127 


-0.309 ± 0.232 


2.21 a 


1.999 


1.663 ±0.134 


1.118 


high-z total 


2.50 


54 


-0.764 ±0.172 


-0.363 ± 0.386 


4.45cr 


2.837 


2.400 ±0.224 


2.099 


Total 


1.67 


128 


-0.457 ±0.102 


-0.332 ±0.211 


4.48 cr 


2.389 


1.974 ±0.127 


1.558 
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Figure 6. Aa/o versus absorption redshift for the previous low-z (open circles), previous high-z (open triangles) and new (open squares) samples. The 
upper panel shows the raw results and lcr error bars. The middle panel shows an arbitrary binning of each sample and the lower panel combines all three 
samples. The redshifts of the binned points are taken as the mean absorption redshift within each bin and the value of Act /a is the weighted mean with its 
associated la error bar. The upper scale is the look-back time to the absorber as a fraction of the age of the universe (Hq = 70 km s -1 Mpc -1 , f! m = 0.3, 
n A = 0.7 -> to = 13.47 Gyr). 



tal number by the number of bins). We plot the weighted mean for 
each bin with the associated 1 a error bar. The lower panel com- 
bines all three samples. Note the expanded vertical scales on the 
lower two panels. 

4.2 Extra scatter at high-z 

In Table 4 we note that the scatter in the total low-z sample is con- 
sistent with that expected from the size of the error bars [i.e. rms « 
(6(Aa/a)) and xt & !]• However, at higher z, Fig. 6 shows sev- 
eral values of Aa/a which are significantly different from one an- 
other, particularly around z a b s ~ 1.9. Indeed, Table 4 shows that 
rms > (5(Aa/a)) and xi > 1 f° r the total high-z sample and 
that the previous high-z sample dominates these statistics. There- 
fore, the weighted mean will not indicate the true significance of 
Aa/ a at high-z. Note that the binned points in the lower panels of 
Fig. 6 will also appear too significant. 

There are several reasons why we expect extra scatter at high 



z. Consider fitting two transitions, arising from different species, 
with significantly different linestrengths (e.g. Al II A 1670 and Ni II 
A 1709 in Fig. 4). Weak components in the high optical depth por- 
tions of the strong transition's profile are not necessary to obtain a 
good fit to the data. Even though the VPFIT \ 2 minimization en- 
sures that constraints on Aa / a derive primarily from the optically 
thin velocity components, these weak components missing from 
the fit will cause small line shifts. The resulting shift in Aa/a is 
random from component to component and from system to system: 
the effect of missing components will be to increase the random 
scatter in the individual Aa/a values. We expect this effect to be 
greater in the high-z sample for several reasons: 

(i) At high z, a larger number of different species are gener- 
ally available for fitting compared to lower z. The range in optical 
depths for corresponding velocity components is therefore signifi- 
cantly larger than for the low-z Mg/Fe II systems (compare Figs. 3 
and 4), complicating determination of the velocity structure (see 
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Figure 7. Reducing the extra scatter at high-z. The upper panel identifies 
the 'high-contrast' systems, i.e. where any of the strong Alii, Sin or Fell 
transitions and any of the weak Cm, Nill or Znll transitions are fitted. 
These systems are expected to exhibit additional scatter. In the middle panel 
the high-contrast systems are removed and the binned low-contrast values 
(solid circles) are compared with the raw values (dotted open circles). In 
the lower panel the high-contrast error bars have been increased by adding 
2.09 X 10 in quadrature to the raw errors to match the observed addi- 
tional scatter. The solid circles are binned values which include the (unal- 
tered) low-contrast systems. 



below for details). Abundance and ionization variations will also 
be more important due to the diversity of species at high-z. 

(ii) The high-z sample is dominated by DLAs. If these have a 
more complex velocity structure, i.e. the number of absorbing com- 
ponents per kms -1 is higher than for the Me/Fell systems, Voigt 
profile decomposition would be more difficult, increasing the scat- 
ter in Aa/a. 

(iii) In the high-z systems we fit more heavy, i.e. low 6, species. 
Absorption features are therefore closer to the resolution of the in- 
strument, increasing the systematic bias against finding the weaker 
components. 

If these effects cause the extra scatter we observe, the most 
affected systems will be those with a large range of transition 
linestrengths. We form such a 'high-contrast' sample from the high- 
z (i.e. z a bs > 1.8) sample by selecting those systems in which any 
of the strong Al II, Si II or Fe II transitions and any of the weak 
Cr II, Ni II or Zn II transitions are fitted. These high-contrast sys- 
tems are marked with a '*' in Table 3. In the top panel of Fig. 7 
these systems (open circles) are delineated from the rest of the 
high-z sample, i.e. the 'low-contrast' systems (solid circles). One 
immediately notes that the high-contrast sample dominates the ex- 
tra scatter. The middle panel in Fig. 7 shows the effect of removing 
the high-contrast systems from the analysis. We give relevant statis- 
tics in Table 4: the high-contrast sample comprises 22 systems and, 
upon removal, xt at high-z drops from 2.179 to 1.207, indicating 
that the extra scatter is indeed due to the effects mentioned above. 

To obtain a more robust estimate of the weighted mean Aa /a 
at high-z, we may increase the individual 1 a errors on Aa/a un- 
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Redshift 

Figure 8. The fiducial sample. The upper panel shows the raw low-z results 
(solid triangles), the raw values for the low-contrast sample (open squares) 
and the high-contrast sample with increased error bars (open circles). The 
middle panel shows a binning of each sample taken separately while the 
lower panel bins all samples together. The weighted mean for this fiducial 
sample is Aa/a = (—0.543 ± 0.116) X 10 -6 and is our most robust 
estimate of Aa/a from the QSO data. 



til Xv = 1 about the weighted mean for the high-contrast sample. 
We achieve this by adding 2.09 x 10~ 5 in quadrature to the er- 
ror bars of the 22 relevant systems. The resulting binned values of 
Aa /a for the entire high-z sample are plotted in the lower panel of 
Fig. 7 together with the raw results (shifted for clarity). The statis- 
tics in Table 4 indicate that the best estimate of Aa/a at high-z is 
Aa/a = (-0.563 ± 0.291) x 10" 5 . Table 4 also gives the statis- 
tics for two other tests: (i) increasing the error bars for all high-z 
systems until \ v — 1 an d (ii) removing points which deviate sig- 
nificantly from the weighted mean (i.e. cr-clipping) in the high-z 
sample until \ v = 1. All tests yield similar results. 

Throughout the rest of this paper, 'fiducial sample' will refer 
to the raw sample (i.e. the values and 1 a errors in Table 3) but 
with increased errors on the high-contrast systems at high-z, as de- 
scribed above. We plot the fiducial sample in Fig. 8 for comparison 
with the raw sample in Fig. 6. Table 4 shows the statistics for the 
fiducial sample: 



Aa/a = (-0.543 ± 0.116) x 10" 



(9) 



is our most reliable estimate of the weighted mean Aa/a over the 
redshift range 0.2 < z a b s < 3.7. This represents 4.7 a statistical 
evidence for a varying a. 



4.3 Temporal variations in a? 

In Fig. 9 we overlay the binned fiducial data from Fig. 8 with sev- 
eral fits to the unbinned values of Aa / a versus z a b s and cosmo- 
logical time. Note that for the constant a/ a = {da /At) /a and 
constant {Aa/Az)/a models we fix Aa/a to zero at z abs = 0. 
The values of \ 2 f° r tne constant Aa/a and constant a/ a fits are 
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X 2 = 135.22 and 134.54 respectively. This indicates that a linear 
increase in a with time, 



a/a = (6.40 ± 1.35) x 10 



-16 -1 

yr 



(10) 



is preferred over a constant Aa/a. We estimate the significance of 
this preference using a bootstrap technique. Each value of Aa/a 
is randomly assigned to one of the measured absorption redshifts 
and the best fitting constant a/ a is determined. Constructing many 
such bootstrap samples yields a distribution of x 2 for the constant 
a/ a fits. In this procedure we construct low- and high-z bootstrap 
samples separately since the distribution of Aa/a in the real low- 
and high-z samples may differ (Section 4.2). The \ 2 distribution 
indicates that 31 per cent of x 2 values were < 135.22 and 15 per 
cent were < 134.54. Thus, the constant a/ a model is preferred 
over the constant Aa/a model only at the 50 per cent confidence 
level. The latter is preferred over a linear evolution of a with z a bs 
at the 63 per cent level. These confidence intervals are only mod- 
estly affected if we analyse the raw sample rather than the fiducial 
sample. Not distinguishing the low- and high-z samples gives lower 
confidence levels, i.e. ^38 per cent for both cases. 

Furthermore, the constant a/ a and (da/dz)/a models 
should be treated with caution. Fixing Aa/a to zero (i.e. fixing a 
to the laboratory value) at z a b s = may not be valid because of the 
potential for spatial variations in a. For example, if a varies over 
~100Mpc spatial scales and our Galaxy resides in a region with a 
slightly larger a than the 'universal mean', then we should expect 
to measure a constant negative Aa/a rather than an evolution of 
a with time. Bekenstein (1979) and Barrow & O'Toole (2001) also 
note that it is difficult to compare values of Aa/a in areas of dif- 
ferent gravitational potential without a detailed theory giving both 
time and space variations of a. If we allow the z a b s = value 
of Aa/a to vary, the linear fit against cosmological time is barely 
altered (see Fig. 9), though somewhat poorly constrained. 



4.4 Angular variations in a? 

The fiducial sample can also be used to search for spatial varia- 
tions in a. If very large scale (~10Gpc) spatial variations exist, 
one might expect Aa/a to be different in different directions on the 
sky. For the fiducial sample, we plot the distribution of Aa/a over 
Galactic coordinates (I, b) in Fig. 10. All redshift information is re- 
moved by taking all absorption systems along a single QSO sight- 
line and collapsing them to a single weighted mean value of Aa /a. 
We also combine independent values from the same QSO in differ- 
ent observational samples (see Section 4.7). The grey-scale indi- 
cates the weighted mean Aa / a and the size of each point scales 
with the significance with respect to the overall weighted mean of 
the fiducial sample, S-0.543 (equation 9). 

We used a x minimization algorithm to find the best-fit dipole 
in the angular distribution of Aa/a. We first constructed a grid of 
directions in equatorial coordinates, (RA, DEC). For each direc- 
tion, a cosine fit to Aa/a as a function of the angular distance, 
<f>, to each absorption system, yields an amplitude for the dipole, 
(Aa/a)d, defined by 



Aa.,, /Aa\ / Aa\ 

W = ( ) + cos< 

a \ a / w V a J d 



(11) 



Fig. 11 shows the best-fit dipole which has an amplitude 
(Aa/a) d = (0.43 ± 0.25) x 10 _s . The pole is in the di- 
rection P A «(RA,DEC) = (5.0hr, -48°) or P A a(/,&) = 
(340°, -3.9°) as marked on Fig. 10. We also mark the anti-pole, 
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Figure 9. Temporal variation in a. (a) The points are the binned values of 
Aa/a for the fiducial sample plotted versus redshift. The dashed line is 
the weighted mean and the solid line is a fit to the unbinned data fixed to 
Aa/a = at z abs = 0. The dot-dashed line is a fit with the z aba = 
value allowed to vary. The values of \ 2 indicate that the constant Aa/a 
model is preferred, (b) Same as (a) but versus fractional look-back time. 
The values of \ 2 indicate that a linear increase in a with time is preferred. 



^Aa, and the cosmic microwave background (CMB) pole and anti- 
pole directions (Pcmb and Acmb, Lineweaver et al. 1996) for 
comparison. 

The limited range of angular separations in Fig. 1 1 severely 
limits the dipole interpretation. Although the dipole amplitude is 
significant at the 1.7 a level, a bootstrap analysis similar to that in 
the previous section indicates that the dipole model is not signif- 
icantly preferred over the constant Aa/a (i.e. monopole) model. 
Bootstrap samples were formed by randomizing the values of 
Aa/a over the QSO sight-line directions. Again, we treated the 
low- and high-z samples separately, combining them only to find 
the best-fit dipole directions. The resulting probability distribution 
for (Aa/a)d indicates that values at >1.7 a significance occur 
40 per cent of the time by chance alone. Analysing the raw sample 
gives a similar confidence level. Therefore, the data do not support 
significant angular variations in a. There is also no evidence for an- 
gular variations in the low-z and high-z samples taken separately. 



4.5 Spatial correlations in a? 

If spatial variations in a do exist at the Aa/a ~ 10~ 5 level then 
one expects to see additional scatter in the raw values of Aa/a 
in Fig. 6. However, as noted previously, the statistics in Table 4 
indicate that the scatter at low z is consistent with the size of the 
error bars and so we have no evidence for spatial variations in a for 
z a b s < 1.8. The additional scatter at higher z may provide evidence 
for spatial variations but this seems unlikely given the discussion in 
Section 4.2. 

We performed a more thorough search for spatial variations 
by calculating the two-point correlation function for a. Consider 
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Figure 10. The distribution of Aa/a in Galactic coordinates (I, b). All fiducial values along a single QSO sight-line (i.e. for many absorption clouds and/or 
repeated, independent QSO observations) are collapsed to a single weighted mean. The value of Aa/a is given by the grey-scale and the size of each point 
scales with the significance, S-o.543, of the departure from the overall weighted mean of the fiducial sample (equation 9). The best-fit dipole has an amplitude 
(Aa/a)<| = (0.43 ±0.25) X 10~ 5 with its pole and anti-pole in the directions P& a and A^ a . However, a bootstrap analysis indicates that the dipole model 
is preferred over the constant Aa/a model at only the 60 per cent confidence level. We plot the CMB pole and anti-pole for comparison. (The on-line version 
of this figure is colour-coded and is also available from http://www.ast.cam.ac.uk/~mim/pub.html) 



two absorption systems with fiducial values of a, ai and ai- The 
two-point correlation function, C12, is then 

(ai - a m ) (a 2 - a m ) 
L-12 = 5 , Ua) 

where a m is the mean value of a over the entire sample. 

In Fig. 12 we plot the correlation function versus the comov- 
ing separation between the absorption systems (e.g. Liske 2003), 
A\- We used a Monte Carlo technique to estimate the mean two- 
point correlation function, C(Ax), and 68 per cent confidence in- 
terval for each bin. Synthetic values of a were drawn randomly 
from Gaussian distributions centred on the measured values with 
1 a widths equal to the fiducial 1 a errors. Where independent ob- 
servations of an absorption system exist (Section 4.7), we centred 
the Gaussian on the weighted mean value and used the 1 a error in 
the weighted mean. For each bin, the mean C( Ax) was found from 
all pairs of absorption systems with comoving separations within 
the bin. In Fig. 12 we plot the mean and rms of this quantity, rep- 



resented by a point and error bar. Fig. 12 shows no evidence for 
significant spatial correlations in a over 0.2-13 Gpc (comoving) 
scales. 

4.6 Comparison with previous analysis 

The majority of velocity structures fitted to the previous samples 
are unchanged from M0 la. For the low-z sample, the raw values of 
Aa/a are similar to those found in MOla, as shown in Fig. 13. In 
MOla, the weighted mean was Aa/a = (-0.70±0.23) x 10~ 5 for 
this sample. The primary reason for the slightly less negative value 
of Aa/a we find here [from Table 4, Aa/a = (-0.51 ± 0.22) x 
10~ 5 ] is the increased q coefficients for the low-z Fell transitions 
(i.e. those with Ao > 2000A in Table 2). For the previous high- 
z sample, Fig. 13 shows that the raw values of Aa/a found here 
are, in individual cases, somewhat different to those found in MOla. 
This is primarily due to the large change in the Fe II A 1608 q coeffi- 
cient (see Section 2.2). However, since the high- z sample is charac- 
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Figure 11. Best-fit dipole to the angular distribution of Aa/a in Fig. 10. 
The pole lies in the direction PA a (l,b) = (340°, -3.9°). The upper panel 
shows the distribution of Aa/a with angular separation from the anti-pole, 
A-Aa - The lower panel shows a binning of the results. The solid line is the 
best-fit cosine with an amplitude (Aa/a) d = (0.43 ± 0.25) X 10~ 5 . 
However, the limited angular coverage of QSO sight-lines clearly under- 
mines confidence in this dipole interpretation. 
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Figure 12. The two-point correlation function for a, C(A\), as a func- 
tion of comoving separation between absorption clouds, A\ (Hq = 
70kms _1 Mpc- 1 , fi m = 0.3, U A = 0.7). The values and errors 
were obtained with a Monte Carlo technique. However, since each value 
of Aa/a may contribute to many bins, the bins are not independent. 
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Figure 13. Comparison with previous analysis in M0 la. The upper panel 
compares the raw values of Aa/a found here with those in M0 la 
[A(Aa/a) = Aa/a — (Aa/a)Mula] f° r me previous low-z (open cir- 
cles) and previous high-z (open triangles) samples. The lower panel shows 
a measure of the significance of this difference: A (Aa/a) normalized by 
the 1 a error in Table 3. For the low-z data we find a slightly more positive 
Aa/a on average because the ^-coefficients for the Ao > 2000A Fell 
transitions have been revised to slightly higher values. The larger scatter 
seen in A( Aa/a) for the high-z systems is due mainly to the large change 
in the Fe II A1608 q coefficient (see Section 2.2). 



terized by many different transitions with a diverse range of q coef- 
ficients, the overall weighted mean Aa /a found in MOla compares 
well with the value found here: Aa/a = (-0.67 ± 0.33) x 10~ 5 
and (-0.67 ± 0.24) x 10~ 5 respectively. 



4.7 Repeated observations 

Several QSOs in our sample have been observed independently by 
two different groups: the new sample contains 7 absorption systems 
which are also (independently) contained in the previous samples 
and two independent spectra of the z^bs = 2.6253 system towards 
GB 1759+7539 are contained in the previous high-z sample. We 
compare the values of Aa/a from the raw sample for these re- 
peated systems in Fig. 14. The lower panel may indicate some ev- 
idence for systematic errors in individual points (i.e. some addi- 
tional scatter may exist). However, the value of \ 2 with respect to 
zero difference is \ 2 = 13.33 which, for 7 degrees of freedom, has 
a 6.4 per cent probability of being exceeded by chance. Evidence 
for additional scatter is therefore rather marginal. 



5 SYSTEMATIC EFFECTS 

The above results represent strong statistical evidence for a smaller 
a in the QSO absorption clouds. The fact that three different sam- 
ples show the same effect leaves little doubt that we do detect small 
line shifts in the data. The central question is now whether these 
line shifts are due to systematic errors or really due to varying a. 
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Figure 14. Repeated absorption systems. The upper panel shows the raw 
values of Aa/a for those systems which have been observed indepen- 
dently by different groups. The absolute difference is shown in the lower 
panel. For 7 degrees of freedom, the indicated value of \ 2 will be exceeded 
6.4 per cent of the time by chance alone, providing only marginal evidence 
for significant additional random errors. 



In MO lb we considered potential systematic effects due to labo- 
ratory wavelength errors, wavelength miscalibration, atmospheric 
dispersion effects, unidentified blending transitions, isotopic and/or 
hyperfine structure effects, intrinsic instrumental profile variations, 
spectrograph temperature variations, heliocentric velocity correc- 
tions, kinematic effects and large scale magnetic fields. Below we 
extend this initial search to include the entire Keck/HIRES sample. 



5.1 Kinematic effects 

In order to measure Aa/a with the MM method one must assume 
that corresponding velocity components of different ionic species 
have the same redshift. Velocity segregation could arise through 
chemical abundance gradients combined with differential velocity 
fields. These effects may also generate departures from Voigt pro- 
files. The low-z Mg/Fe II systems would be more sensitive to any 
velocity segregation between species due to the simple arrangement 
of q coefficients (see Fig. 1). For the high-z systems, the effect on 
Aa/a for an individual system would be minimized due to the di- 
versity and complicated arrangement of q coefficients. 

Inhomogeneous photoionization throughout the absorption 
complex is also a possibility. With the exception of Mg I and Al III, 
the species in our analysis are singly ionized, have very similar 
ionization potentials, IP + , all above the Lyman limit and the cor- 
responding neutral ions have similar ionization potentials, IP - , all 
below the Lyman limit (see Table 2). Thus, the absorption from dif- 
ferent ionic species, associated with a given velocity component, 
should arise co-spatially. However, variations in the incident ra- 
diation field could cause weak components in the singly ionized 
species to appear very strong in Mg I or Al III. Spurious shifts in 
Aa/a could be introduced if a strong Mgl or Aim component 



was mistakenly tied to a nearby (|At>| ~ 0.5 kms~ ) singly ion- 
ized component 7 . 

However, as we pointed out in MO lb, these kinematic shifts 
in Aa/a are themselves random over a large ensemble of random 
sight-lines. If kinematic effects are significant, they act only to in- 
crease the scatter in the values of Aa/a beyond what is expected 
on the basis of the statistical errors alone. The value of xl for the 
low-z sample in Table 4 is «1 and so the low- z data present no evi- 
dence for significant kinematic effects. One may interpret the extra 
scatter in Aa/a observed at high-z as evidence for such effects, 
though we discuss other, more likely, mechanisms for producing 
this scatter in Section 4.2. 

Despite the random nature of kinematic effects over a large 
number of velocity components and absorption systems, we discuss 
below the contributory factors to potential effects on small velocity 
scales for individual velocity components. 



5.1.1 Small-scale velocity structure and equilibrium in QSO 
absorbers 

Detailed kinematic studies of QSO absorbers indicate that a 
disk+halo model provides a reasonable description of the absorp- 
tion systems (Briggs et al. 1985; Lanzetta & Bowen 1992; Wolfe & 
Prochaska 2000b; Churchill & Vogt 2001). Alternative models ex- 
ist, including multiple merging clumps bound to dark matter halos 
(Haehnelt, Steinmetz & Rauch 1998; McDonald & Miralda-Escude 
1999; Mailer et al. 1999) and outflows from supernovae winds 
(Nulsen, Barcons & Fabian 1998; Schaye 2001). For disk+halo 
models of Mg/Fe II systems, the disk component is strongly sat- 
urated and spread over small velocity scales whereas the halo com- 
ponent is broadly spread in velocity space and causes the lower col- 
umn density absorption. The constraints on Aa/a in the Mg/Fe II 
systems are therefore dominated by velocity components arising 
in the outer parts of galaxy halos. The constraints on Aa/a from 
DLAs arise either from low abundance, hence unsaturated, species 
(e.g. Ni II, Cr II, Zn II), or from optically thin components flanking 
the saturated ones (e.g. Fell A1608 or Alii A1670). 

However, the large-scale properties of the absorbing gas have 
no influence on estimates of Aa/a. Galactic rotation or large-scale 
galactic winds are unimportant. The empirical fact that we find 
excellent agreement between the redshifts of individual velocity 
components down to < 0.3 km s -1 illustrates this. Any contribu- 
tion to the scatter in Aa/a comes only from the detailed proper- 
ties of gas on velocity scales typical of the observed b parameters, 
i.e. < 5kms _1 . 

Given the importance of small-scale properties in the deter- 
mination of Aa/a, it is relevant to ask whether the gas is in dy- 
namic equilibrium on these scales. For the low-z Mg/Fe II systems, 
a recent detailed study (Churchill, in preparation) suggests that the 
column density ratio 7V(Mg Il)A/V(Fe II) appears not to change sys- 
tematically across an absorption complex. A similar chemical uni- 
formity is found for the majority of high-z systems (Prochaska 
2003). This implies that, if not in equilibrium, the gas is well 
mixed, i.e. absorption lines of different species arise co-spatially. 



7 We find no additional scatter in Aa/a for systems containing Mg I or 
Aim and so have no evidence for such effects. We discussed the specific 
case of Aim in M0 lb, finding that its velocity structure followed very 
closely that of the singly ionized species. See also discussion in Wolfe & 
Prochaska (2000a) 
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Photoionization models also suggest gas in photoionization equi- 
librium with an ambient extra-galactic UV background. If this is 
correct, local equilibrium may be valid. Any redshift evolution of 
the number of absorption lines per unit redshift interval, over and 
above that expected due to cosmology alone, could be explained by 
cosmological evolution of the integrated background UV flux (e.g. 
Rao & Turnshek 2000). 

Thus, there is evidence to suggest that no gross changes in 
physical conditions occur over large velocity scales across an ab- 
sorption complex. This may imply that we should not expect to find 
abundance variations and non-equilibrium on small scales. If this is 
incorrect, one expects departures from a Maxwellian velocity dis- 
tribution in the absorbers, the assumption of which is inherent in 
fitting Voigt profiles. We note in passing that no evidence currently 
exists for non- Voigt profiles (e.g. Outram et al. 1999), despite ex- 
pectations from hydrodynamic simulations of large-scale structure 
formation (Outram et al. 2000). 



5.1.2 Comparison with the ISM 

In contrast to the previous section, analogy with the interstellar 
medium (ISM) in our own Galaxy suggests that non-equilibrium 
could apply on very small scales. Andrews, Meyer & Lauroesch 
(2001) used stars in a background globular cluster, M92, to probe 
the kinematics on scales defined by the separation between the 
lines of sight at the absorber. They find significant variations in 
Na I column densities in the ISM on scales as small as 1600 AU (or 
~ 0.01 pc). Even smaller-scale details come from measurements of 
temporal variation of Na I and K I absorption lines, implying non- 
equilibrium scales ~10-100 AU (Crawford et al. 2000; Lauroesch, 
Meyer & Blades 2000; Price, Crawford & Barlow 2000). Although 
these ISM cloud sizes are small compared to estimates for individ- 
ual QSO absorption cloud components, ~ 10-100 pc (Churchill & 
Vogt 2001), the characteristic size for a QSO continuum emission 
region may be ~ 10~ 3 pc. The lines of sight to opposite edges of 
the QSO therefore probe similar size scales as the Galactic ISM 
studies. However, it should be noted that the gas densities are very 
different and so the comparison should be treated with caution. 

The ISM analogy may also allow us to estimate an upper limit 
on the velocity segregation in QSO absorption clouds. Some the- 
ories argue that gravity may be important for cloud confinement 
on small size scales in the Galactic ISM (e.g. Walker & Wardle 
1998). If gravity plays a significant role in QSO absorption sys- 
tems on similar scales, we could apply a simple stability condition 
on the velocity dispersion, v 2 = GM/R, where M and R are the 
cloud mass and radius. Estimates for individual cloud sizes vary but 
adopting R = 10 pc and M = 30 M (Churchill & Vogt 2001) we 
find \v\ < 0.1 kms -1 . This provides an upper limit on the veloc- 
ity shift between different species. For one single Mg/Fe II veloc- 
ity component, this translates into an error on an individual Aa/a 
measurement of roughly |Aa/a| ~ 0.3 x 10~ 5 . However, this 
would be randomized over ~100 observations (and over many ve- 
locity components), producing a maximum overall effect which is 
20 times smaller than that observed. 



5.2 Line blending 

We distinguished between random and systematic line blending in 
MOlb. We discuss both cases and their effect on Aa/a below. 



5.2.1 Random blends 

Random blends can occur when many absorption clouds lie along 
the line of sight to the QSO, including any interstellar material and 
the Earth's atmosphere. In general, random blends can only have a 
random effect on Aa/a and so can not have caused the systemati- 
cally non-zero values observed. Two distinct categories of random 
blend can be identified: 

(i) Strong blends. Lyman-a forest absorption lines would cause 
numerous strong blends. However, if MM transitions fell into the 
Lyman-a forest region, we generally did not fit them. If the veloc- 
ity structure of one of our MM transitions is obviously affected by 
a random blend then we modify our fit for that transition according 
to the identity of the interloper. If the interloper is a MM transi- 
tion (i.e. one with a precise laboratory wavelength listed in Table 
2), and its velocity structure can be constrained using other asso- 
ciated MM transitions at the interloper's redshift, then we include 
the interloper in the fit. In general, the laboratory wavelengths of 
non-MM transitions (e.g. Civ AA1548 and 1550) are not known 
to high precision. Therefore, if the interloper and/or the associated 
transitions are not MM transitions, possible errors in the labora- 
tory wavelengths could have a significant effect on Aa/a. In these 
cases, and if the velocity structure of the interloper was both simple 
and clear, we freed its redshift parameters, not tieing it to any asso- 
ciated transitions. In all other cases, we either masked the blended 
data out of the fit or simply rejected that MM transition altogether. 

(ii) Weak blends. It is possible that weak, random interlop- 
ers exist in our data but were not identified in our analysis. For 
example, numerous broad (FWHM = 20-2000 kms -1 ), weak 
(equivalent width < 0.1 A), diffuse interstellar absorption bands 
have been discovered in stellar spectra (e.g. Herbig 1975; Jen- 
niskens & Desert 1994). The narrower of these lines may blend 
with some velocity components of individual MM transitions, caus- 
ing slight, apparent line shifts. However, since our absorption sys- 
tems lie at a range of redshifts and since many different transitions 
are fitted to arrive at a final value of Aa/a, it is unlikely that such 
weak interlopers have significantly affected our results. 



J. 2. 2 Systematic blends 

A systematic blend occurs when two ionic species are in the same 
absorption cloud and have transitions with similar rest wavelengths. 
Such a blend could mimic a systematic shift in a. We have explored 
two approaches to this problem: 

(i) In MOlb we attempted to identify candidate blends by 
searching atomic line databases for transitions lying close to the 
laboratory wavelengths of the MM transitions. Simulations of var- 
ious blends indicated that, in order to cause a significant and sys- 
tematic shift in a, the interloping transition must lie within \Av\ w 
lOkrns -1 of the 'host' line. No possible blending transitions sat- 
isfying this criterion could be found in Moore (1971) or the Vienna 
Atomic Line Database 8 (VALD, Piskunov et al. 1995; Kupka et al. 
1999). 

We have now extended our search criterion to include potential 
interlopers lying up to 100 kms -1 away from the host MM transi- 
tions. We have identified two possible interlopers lying to the red 
of the Zn II A2026 transition: 



Available at http://www.astro.univie.ac.at/~vald 
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(a) Cm A2026.269 (Av w 19kms _1 ) has an oscillator 
strength / = 0.0047 (VALD; correcting a typographical er- 
ror in MOlb). Using the largest value of N(Cr II)/ AT (Zn II) from 
our data, simulations of this blend suggest a negligible effect on 
Aa/a (see MOlb for details). 

(b) Mgl A2026.477 (At; « 50km s _1 ) has an oscillator 
strength / = 0.1154 (Morton 1991) and was apparent in six of 
the high- z systems. However, in four cases the velocity structure 
is narrow enough that the Mg I interloper does not overlap and 
blend with the Znll line. In the remaining two cases we fitted 
the Mg I interloper and left the redshift parameters free (i.e. we 
did not tie the Mg I velocity structure to that of the other MM 
transitions). This has the effect of removing any constraints the 
blended Zn II components may have on Aa/a. We checked the 
resulting values of Aa/a against those found after completely 
masking out the blended portion of the Zn II velocity structure, 
finding them to be consistent. Also, completely removing the 
Znll A2026 transition from the analysis of these systems pro- 
duced very little change in Aa/a. 

Note that the atomic line databases searched are not complete 
and that we have not considered molecular interlopers. We there- 
fore compliment the treatment above with the following test. 

(ii) To asses the significance of any unknown blends, we may 
separately remove most transitions from the analysis and re-fit to 
find a new value of Aa/a. If the new and old values of Aa/a 
are significantly different then one interpretation is that the transi- 
tion removed is affected by an unknown systematic blend. In MOlb 
we found no evidence for blends in this way. We update this line- 
removal analysis in Section 5.7.1 and reach the same conclusion. 

In summary, we find that neither random or systematic line 
blending has significantly affected the values of Aa/a. 

5.3 Wavelength miscalibration 

The QSO CCD images are wavelength calibrated by comparison 
with ThAr emission lamp spectra which are generally taken both 
before and after each QSO exposure. A polynomial wavelength so- 
lution is obtained by centroiding a standard set of ThAr lines which 
are strong and (appear to be) unblended in the particular ThAr lamp 
in use. If one or more of the ThAr lines in this standard set are 
misidentified or, in a more subtle way, lead to long range miscal- 
ibrations of the wavelength scale, then a systematically non-zero 
Aa/a would be found. This effect would be particularly prob- 
lematic in the low-z absorption systems where all values of Aa/a 
would be affected in the same way by such a low-order distortion 
of the wavelength scale (see Fig. 1). 

5.3.1 The ThAr calibration test 

The effect of any ThAr miscalibrations on the values of Aa/a can 
be estimated directly using a technique proposed and demonstrated 
in MOlb: we analyse sets of ThAr emission lines in the same way 
we analyse sets of QSO absorption lines. If no calibration error ex- 
ists then one should obtain (Aa/a)ThAr = for each absorption 
system. 

We provide a detailed description of this ThAr test in MOlb, 
the main steps being as follows. First we select several independent 
sets of ThAr emission lines from sections of spectra which corre- 
spond to the observed wavelengths of the QSO absorption lines. 
Each set of lines therefore corresponds to, and has been selected in 



an analogous way to, one QSO absorption system. Selecting several 
(typically 3-7) different sets of ThAr lines allows a more robust es- 
timate of (Aa/a) T hAr for each absorption system and provides a 
cross-check on its 1 a uncertainty. 

We fit Gaussian profiles to the emission lines with the same 
modified version of VPFIT used to fit Voigt profiles to the QSO 
absorption lines. In order to relate any measured ThAr line shifts 
directly to a value of (Aa/a) T hAr, we assign to each ThAr line 
the q coefficient corresponding to the QSO absorption line falling 
in that part of the spectrum. That is, we treat the ThAr lines as if 
they were QSO absorption lines and derive a value of (Aa/a)ThAr 
using the MM method. 



5.3.2 Results and discussion 

ThAr spectra were readily available for 118 of the 128 absorption 
systems in our total sample. We applied the above ThAr test to these 
spectra and obtain the results presented in Fig. 15. Comparing the 
QSO results (raw sample) in panel (a) with the ThAr results in panel 
(b) we immediately see that miscalibrations of the wavelength scale 
did not systematically affect the overall value of Aa/a. This is 
particularly emphasized in panel (e) which compares binned raw 
values of Aa/a and (Aa/o^ThAr- We tabulate various statistics 
for the ThAr test in Table 4 for direct comparison with the QSO 
results. Notably, the overall weighted mean value, (Aa/a)xhAr = 
(0.4 ± 0.8) x 10 -7 , is two orders of magnitude below the observed 
value of Aa/a. 

Panel (c) of Fig. 15 shows an expanded view of the ThAr re- 
sults which clearly demonstrate a similar 'extra scatter' to the high- 
z QSO sample (see Section 4.2). The extra scatter has a similar ori- 
gin as for the high-z QSO results: close inspection of any portion 
of ThAr spectrum reveals a plethora of weak emission lines and 
so the fit to any given strong ThAr emission line will be affected 
by weak blends. The ThAr lines are of such high S/N that signif- 
icant extra scatter ensues. This effect would produce uncorrelated 
deviations in (Aa/a)ThAr from point to point, as observed. An- 
other interpretation may be that the extra scatter presents evidence 
for ThAr line misidentifications. However, if these were important, 
one would expect a correlation between (Aa/a) T hAr and Aa/a. 
We observe no such correlation. 

The ThAr results in Fig. 15 also illuminate discussion of 
potential instrumental profile (IP) variations and their effect on 
Aa/a. Valenti, Butler & Marcy (1995) have derived the HIRES IP 
for several positions along a single echelle order, finding that the IP 
is slightly asymmetric and that this asymmetry varies slightly along 
the order. However, to a first approximation, the ThAr and QSO 
light paths through the spectrograph are the same. The ThAr re- 
sults in Fig. 15 therefore strongly suggest that IP asymmetries (and 
variations thereof) do not contribute significantly or systematically 
to Aa/a. One possibility is that the ThAr and QSO light follow 
very different and possibly wavelength dependent paths through the 
spectrograph to the CCD. We explore this possibility in Section 5.6 
in the context of atmospheric dispersion effects. 

We conclude this section by emphasizing the reliability of 
the wavelength scale as derived from the ThAr spectra: we di- 
rectly measure the effect of any ThAr miscalibrations on Aa/a 
and find it to be negligible. One objection to the reliability of the 
MM method has been that fitting transitions which fall on different 
echelle orders may lead to systematic errors (Varshalovich et al. 
2000; Ivanchik et al. 2002). The ThAr test in Fig. 15 (and figure 2 
in MOlb) clearly demonstrates that these concerns are unfounded. 
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Table 5. The percentage terrestrial isotopic abundances of the atoms used in 
our analysis (Rosman & Taylor 1998). The second column shows the mass 
number, A, of the isotope with the highest abundance and AA is defined 
relative to this (negative values representing lighter isotopes). 
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Figure 15. Comparison between the raw QSO results [panel (a)] and the 
ThAr calibration error [panel (b)] for the previous low-z (open circles), 
previous high-z (open triangles) and new (open squares) samples. The 
(Aa/a)xtAr scale is expanded by a factor of 10 in panel (c). Panel (d) 
shows binned values for the three sub-samples and panel (e) compares the 
binned values for the total ThAr (solid circles) and raw QSO (open cir- 
cles) samples. The weighted mean of the ThAr results is (Aa/aJxhAr = 
(0.4 ±0.8) X 10~ 7 . 



5.4 Spectrograph temperature variations 

The refractive index of air within the spectrograph depends on tem- 
perature (and also on pressure, but this is a smaller effect). If the 
temperature during the QSO exposure, Tqso, is ~ 15 K higher 
than that during the ThAr calibration exposure, TxhAr, the wave- 
length scale applied to the QSO frame would be distorted, resulting 
in Aa/a ~+lx 1CT 5 for a typical Mg/Fell system (see MOlb 
for details). Note that a systematically non-zero Aa/a can only 
be introduced for the Mg/Fe II systems if Tqso is systematically 
higher or lower than TrhAr- 

Such an effect can be minimized if the ThAr calibration expo- 
sures are taken immediately before and/or after the QSO exposures 
or if the spectrograph is temperature stabilized (e.g. VLT/UVES). 
The former was generally the case for our observations. We have 
used image header information to calculate AT = (Tqso) — 
(T T hAr) for each QSO in each sample, where the average is taken 
over all QSO and ThAr exposures. We find mean values of AT = 
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0.04 ± 0.02 K and AT = 0.2 ± 0.1 K for the previous low- and 
high-z samples, and AT = 0.013 ± 0.006 K for the new sample. 
It is therefore clear that our values of Aa/a are not affected by 
spectrograph temperature variations. 



5.5 Isotopic and hyperfine structure effects 

As stated in Section 2.2, the only species for which full isotopic 
and hyperfine structures are known are Mg I, Mg II and Al III. We 
obtained an estimate of the Si II isotopic structures by scaling the 
Mg II A2796 structure by the mass shift, equation 7. In this section 
we investigate systematic errors that could result from ignorance 
of the isotopic/hyperfine structures for other species and from any 
evolution of the isotopic ratios and/or hyperfine level populations. 



5.5.1 Differential isotopic saturation 

For the transitions of Cr II, Fe II, Ni II and Zn II, the laboratory 
wavelengths in Table 2 are composite values only. The composite 
wavelengths are only strictly applicable in the optically thin regime 
(linear part of the curve of growth). Consider an absorption line 
with several isotopes. As the dominant (highest abundance) isotope 
saturates, the weaker isotopes will have increased influence on the 
fitted line centroid. This could lead to systematic errors in Aa/a. 

Table 5 shows the terrestrial isotopic abundances for the atoms 
used in the MM method. Note that the abundance is quite dis- 
tributed and asymmetric for Ni, Zn and, to a lesser extent, Cr. If 
the isotopic components of the transitions of these species were 
widely spaced in wavelength then differential isotopic saturation 
may be an important effect. However, the Cr, Ni and Zn lines are 
always weak in our QSO spectra, typically absorbing only ~ 20 per 
cent of the continuum. Using only the composite wavelengths for 
transitions of these species is therefore justified. 

Although the Fe abundance is strongly centred on the Fe 
isotope, the Fe II transitions in our QSO spectra are often saturated. 
The possible effect of differential isotopic saturation for Fe should 
therefore be tested. The mass isotopic shift (equation 7) for Fe tran- 
sitions should be > 5 times smaller than for those of Mg. There- 
fore, using only composite wavelengths of Mg and Si transitions in 
our analysis (instead of the full isotopic structures) should place a 
firm upper limit on the possible effect from Fe II. 

We have conducted such a test with our total sample: we use 
composite wavelengths for the Mg I, Mg II and Si II isotopic struc- 
tures and refit our spectra to find values of (Aa/a) com p. We also 
use composite wavelengths for the Al III hyperfine structures (see 
Section 5.5.2). We provide a detailed discussion of the expected 
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Figure 16. Upper limits on the effects of differential isotopic saturation and 
isotopic abundance evolution. The top panel compares the binned values 
of Aa/a derived using only the composite laboratory wavelengths (solid 
triangles) from Table 2 and the values from the raw sample (dotted circles, 
slightly shifted for clarity). The solid circles in the middle panel are derived 
using only the isotopic components of the Mg I, Mg II and Si II transitions 
with the highest terrestrial abundance. The lower panel shows the differ- 
ences between these composite and single isotope values of Aa/a and the 
values from the raw sample. 



magnitude of A(Aa/a) = (Att/a) comp — Aa/a in MOlb and 
find values of (Aa/a) mmp for the previous low- and high-z sam- 
ples. The top panel of Fig. 16 compares values of (Aa/a) comp 
(solid triangles) for the total sample with those found using the full 
isotopic/hyperfine structures (dotted circles). We see that isotopic 
saturation has only a small effect on Aa/a, even for the Mg/Fe II 
systems at low-z. This is confirmed when comparing the statistics 
for this test in Table 4 with the values from the raw sample. It is 
therefore clear that differential isotopic saturation of Fe II transi- 
tions is unlikely to have significantly affected our results. 



J. 5.2 Differential hyperfine saturation 

Hyperfine splitting of energy levels occurs only in species with odd 
proton or neutron numbers (note that most isotopes in Table 5 are 
even). Different hyperfine components will have different transi- 
tion probabilities but the composite wavelength of a line will be 
unchanged by the splitting (i.e. the centre of gravity of the hyper- 
fine components is constant with increased splitting). However, a 
similar differential saturation effect will occur for the hyperfine 
components as for the isotopic components discussed in Section 
5.5.1. 

By far the most prominent hyperfine structure is that of Al III 
which is clearly resolved in the laboratory experiments of Gries- 
mann & Kling (2000) and which we take into account in our fits to 
the QSO data (see Section 2.2). The high-z composite wavelength 
results in Fig. 16 therefore provide an upper limit on the effect of 
hyperfine saturation effects for all other transitions. The only other 
known hyperfine structures are those of the 25 Mg transitions. These 



are so closely separated that they would have a negligible effect on 
Aq / a (especially considering the small effect of differential iso- 
tonic' saturation in Fig. 16). In MOlb we ruled out significant effects 
on Aa/a from saturation effects in other (odd) species due either 
to their low relative abundances, low magnetic moments or the na- 
ture of the ground- and excited-state wavefunctions. All the above 
arguments imply a negligible effect due to differential hyperfine 
saturation. 

J. 5.3 Isotopic abundance variations 

We assumed the terrestrial isotopic abundances in Table 5 when 
fitting Mg and Si absorption lines. However, if the isotopic abun- 
dances in the QSO absorbers are different, small apparent shifts 
in the absorption lines would be introduced, potentially mimicking 
a non-zero Aa/a. This effect will be greatest for changes in the 
isotopic abundances of Mg I, Mg II and Si II since the mass shift 
(equation 7) implies these ions will have the largest isotopic sepa- 
rations. 

Observations of Mg (Gay & Lambert 2000) and theoretical 
estimates for Si (Timmes & Clayton 1996) in stars clearly show 
a decrease in the isotopic abundances with decreasing metallic- 
ity. For example, at relative metal abundances [Fe/H] ~ —1, 
25 Mg/ 24 Mg « 26 Mg/ 24 Mg « 0.1, which is about 30per cent 
below terrestrial values (see Table 5). Theoretical estimates sug- 
gest even larger decreases (Timmes, Woosley & Weaver 1995). The 
metallicities of our absorption systems are all likely to be sub-solar: 
the Mg/Fe II systems (previous low-z sample) have Z — —2.5- 
0.0 (Churchill & Le Brun 1998; Churchill et al. 2000a,b) and the 
DLAs in the previous high-z sample have Z « —1.0 (Prochaska 
& Wolfe 1999, 2000). Therefore, we expect significantly lower iso- 
topic abundances of 25 ' 26 Mg and 29 > 30 Si in our QSO absorption 
systems. 

Considering the above, refitting the QSO spectra using only 
the strong isotopes of Mg and Si (i.e. 24 Mg and 28 Si) should place 
an upper limit on the effect of possible isotopic abundance vari- 
ations on Aa/a. We first demonstrated this test on the previous 
samples in MOlb and have now applied it to the total sample. We 
find similar results. 

Our new results, (Aa/«)i S o (solid circles), are compared with 
the values from the raw sample (dotted circles) in the middle panel 
of Fig. 16. The lower panel also shows the size of the correction, 
A(Aa/a) = (Aa/a)i so — Aa/a. Fig. 16 clearly demonstrates 
that strong isotopic abundance evolution in Mg and Si can not 
explain our results. Indeed, Aa/a becomes more negative when 
we remove this potential systematic effect. Note that removing 
25 ' 26 Mg affects only the low-z Mg/Fe II systems whereas removing 
29 ' 30 Si affects only the high-z points. The largest effect is for the 
low-z Mg/Fe II systems. This is expected because (a) the isotopic 
separations for Mg transitions are larger than those for Si transi- 
tions (equation 7), (b) the relative terrestrial abundances of 25 Mg 
and 26 Mg are larger than 29 Si and 30 Si (Table 5) and (c) a large 
number of different species are fitted in the high-z DLAs so a sys- 
tematic effect in just one species (in this case Si II) is unlikely to 
significantly affect Aa/a. 

The above test estimates the maximum effect of strong iso- 
topic abundance evolution for Mg and Si transitions. Although we 
emphasize that the effect for heavier species should be reduced by 
the mass shift, we point out that isotopic (and hyperfine) structures 
for transitions of other species are not known. It is possible that 
large specific isotopic shifts exist for some transitions. If this is 
combined with strong isotopic ratio evolution then small apparent 
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line shifts would be measured in the QSO spectra. However, this 
possibility is unlikely given the results of removing single transi- 
tions or entire species in Section 5.7.1. 

One possibility is that the isotopic abundances do not fol- 
low the above expected trend based on the Galactic observations 
and theoretical models. For example, higher 25,26 Mg abundances 
in the absorption clouds would mimic Aa/a < in the low-z 
Mg/Fe II systems. The results above allow us to estimate the Mg 
isotopic ratios required to explain the observed Aa/a at low-z. 
The approximate shift in Aa/a will be A(Aa/a) ~ /3Au> s /u s 
where Au> g /u>g is the fractional change in the abundance weighted 
mean transition wavenumber. Applying this relation to the test re- 
sults above, we find that the proportionality constant (3 = 12.2 
using an average A(Aa/a) — —0.445 x 1CP° (from Table 4) 
and Acjg/cjg = —3.6 x 1CT 7 , the fractional difference between 
the 24 MgII and composite Mgll A2796 wavenumbers. To shift 
Aa/a to an average of zero in the low-z sample, we require 
A(Aa/a) = +0.539 x 10~ 5 . Using the (3 just derived, this 
implies Awg/wg = 4.4 x 10~ 7 which could be achieved if the 
( 25 Mg+ 26 Mg)/ 24 Mg abundance ratio is ^0.58 in the absorption 
clouds (compared to the terrestrial value of 0.27). 

5.5.4 Hyperfine level population variations 

Once again, consider the prominent hyperfine structure of the Al III 
transitions. If the populations of the ground-state hyperfine lev- 
els are not equal in the absorption clouds, as they are in the lab- 
oratory measurements of Griesmann & Kling (2000), we would 
measure small apparent shifts for the Aim transitions. In MOlb 
we pointed out that a lower bound on the relative population is 
set by interaction with CMB photons, exp[— /iAo>//cbTcmb(z)], 
where Alj ~ 0.5 cm -1 is the hyperfine splitting for the Aim 
AA1854 and 1862 transitions. At z abs ~ 2.0, the relative popula- 
tion is ~0.9, leading to a shift in the line centroid of ~0.01 cm -1 
(~5x 10~ 4 A). This corresponds to Aa/a = +1 x 10~ 5 for the 
Al III lines alone. Clearly, including other transitions will reduce 
the effect on Aa /a for a particular system. 

Two caveats significantly lessen concern in this case: (i) colli- 
sional excitation processes should drive the relative populations to 
equality; (ii) removing either or both of the Al III transitions from 
our analysis has a very small effect on Aa/a (see Section 5.7.1). 
We therefore consider systematic errors from variations in the hy- 
perfine level populations to be negligible. 

5.6 Atmospheric dispersion effects 

Below we describe how atmospheric dispersion could have dis- 
torted the wavelength scale and instrumental profile (IP) of the 
QSO spectra. This was the largest potential systematic effect iden- 
tified in MOlb for the previous samples. We model the potential 
effects on our QSO spectra in Section 5.6.1, significantly improv- 
ing the model in MOlb. We search for such effects in the QSO data 
in Section 5.6.2, finding that the data themselves suggest negligible 
effects on Aa/a from atmospheric dispersion. 

5. 6. 1 Modelling possible distortions from atmospheric 
dispersion 

The Keck/HIRES was only fitted with an image rotator in 1996 Au- 
gust. Therefore, prior to this time, the spectrograph slit could not 
be held perpendicular to the horizon during observations and so the 



atmosphere will have dispersed QSO light with some component 
across the slit (see figure 5 of MOlb). All of the 27 absorption sys- 
tems in the previous low-z sample were observed before 1996 Au- 
gust. Similarly, 1 1 systems in the previous high-z sample and half 
(39) of the new sample were observed 'pre-rotator' . Therefore, for 
the total sample, 77 of the 128 absorption systems could be affected 
by atmospheric dispersion effects. 

Before the image rotator was installed, the angle of the 
Keck/HIRES slit to the vertical, 9, was directly related to the zenith 
angle by 6 = £. That is, the slit lies with its axis along the horizon 
at £ = 90°. The component of atmospheric dispersion along the 
spectral direction of the slit (i.e. across its axis) therefore increases 
with increasing £ because atmospheric dispersion and 9 increase 
with £. This will lead to two effects on the QSO wavelength scale: 

(i) Compression. Due to the angular separation of different 
wavelengths as they enter the spectrograph, the QSO spectrum will 
be distorted relative to the ThAr calibration frames. Consider two 
wavelengths, A2 > Ai, entering the spectrograph slit separated by 
an angle Atp in the vertical direction. The spectral separation of the 
two wavelengths in the extracted, calibrated spectrum will be 

aAib sin 9 

AA w A 2 - Ai — - , (13) 



where a is the CCD pixel size in angstroms and S is the pro- 
jected slit width in arcseconds per pixel (for HIRES, 8 = 
O^STper pixel). We checked the accuracy of this equation with 
ZEMAX models of the Keck/HIRES provided by S. Vogt. If Atp 
then the spectral separation between any two wavelengths will 
decrease, i.e. the spectrum is compressed 9 . If Atp is due to atmo- 
spheric dispersion then it is a function of the atmospheric condi- 
tions and can be estimated using the refractive index of air at the 
observer and the zenith angle of the QSO (e.g. Filippenko 1982). 
Note that equation 13 assumes that the seeing profiles at Ai and A2 
are not truncated by the slit edges. 

(ii) Wavelength dependent IP asymmetry. If tracking errors or 
seeing effects do cause profile truncation at the slit edges then the 
truncation will be asymmetric and this asymmetry will be wave- 
length dependent. For example, the optical design of Keck/HIRES 
is such that a blue spectral feature will have its red wing truncated 
and a red feature will be truncated towards the blue. The asym- 
metry and the severity of its wavelength dependence will increase 
with increasing Atp and with 9 (and therefore, £)• N° te that when 
we centroid absorption features in the QSO spectra, we will find 
larger wavelength separations due to this effect, i.e. a positive term 
is effectively added to the right-hand-side of equation 13. 

For the low-z Mg/Fe II systems, these two effects will have an 
opposite effect on Aa/a: the compression will produce Aa/a < 
and the wavelength dependent asymmetry will result in Aa/a > 
0. It is important to note that both effects rely mainly on the same 
parameters for each observation (i.e. Atp and 9 — £) and so their 
relative strength is fixed by our model. The one free parameter is the 
seeing, which determines the strength of the wavelength asymmetry 
effect and not the compression effect. We discuss this below. 

For each QSO observed without the image rotator we have 
calculated an effective seeing profile at the spectrograph slit. In 

9 Note that in MOlb we argued that this effect amounted to an expansion 
of the spectrum rather than a compression. The change in sign here is due to 
the optics of the CCD camera on Keck/HIRES which reverses the image on 
the CCD (T. Bida & S. Vogt, private communication). This has now been 
checked with ZEMAX models of the HIRES. 
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general, each QSO spectrum is the combination of several expo- 
sures and each exposure was taken at a different £. For each ex- 
posure we assume typical observing conditions at the telescope 10 
and use the mean value of £ calculated from the recorded image 
headers to obtain the relative angular separations, Aip, between all 
observed wavelengths. Then, to reconstruct the full dispersion pat- 
tern relative to the slit edges, we assumed that light of wavelength 
A = 5500 A was centred on the slit axis. This corresponds to the 
mean wavelength of the response curve for the acquisition camera 
used to guide Keck/HIRES on the QSO. 

For each exposure we assumed a seeing of FWHM = 0'.'75 
to produce a Gaussian intensity profile at the slit for each observed 
wavelength. Although this is a realistic assumption, we have no 
detailed information about the seeing conditions for each obser- 
vation 11 . Tracking errors will broaden the intensity profile so the 
above assumption is effectively a conservative one: if FWHM > 
Ci.'75 then the effect of wavelength dependent IP asymmetry on 
Aa/a will reduce because each wavelength will illuminate the slit 
more uniformly. 

For each QSO observation and at each observed wavelength 
we average the slit intensity profiles for all exposures. For each ob- 
served wavelength (i.e. each observed MM transition) a high S/N, 
high spectral dispersion synthetic spectrum is constructed based on 
the Voigt profile fits to the QSO data. We convolve each spectrum 
with the corresponding intensity profile, truncated by the slit edges, 
and then we convolve with a Gaussian instrumental response of 
width FHWM w 2.2 km s -1 (Valenti et al. 1995). We vary this 
parameter slightly depending on the observed wavelength so as to 
ensure constant velocity resolution to match the real QSO data. Fi- 
nally, we re-sample the simulated spectra to match the real QSO 
spectral dispersion. 

The above procedure provides synthetic QSO spectra with 
Aa/a — but which contain effects due to our model of atmo- 
spheric dispersion. We fit these spectra and determine (Aa/a) a( j, 
an estimate of the effect of atmospheric dispersion on the values of 
Aa/a derived from the real QSO spectra. The results are plotted 
in the top panel of Fig. 17 for the 77 systems observed without the 
image rotator. In the low-z Mg/Fe II systems, atmospheric disper- 
sion causes Aa/a < (i.e. the compression effect dominates the 
wavelength dependent IP asymmetry effect) and so, once the effect 
is modelled and removed, Aa/a increases. At high z not all val- 
ues of Aa/a react the same way to atmospheric dispersion since 
different transitions are fitted in each absorption system. However, 
on average, the high-z values of Aa/a decrease after removing 
atmospheric dispersion effects. 

The middle panel of Fig. 17 shows the corrected values of 
Aa/a, i.e. (Aa/a) adc = Aa/a — (Aa/a) ad (solid squares), 
including those systems not affected by atmospheric dispersion ef- 
fects [(Aa/a) ad = 0, open squares]. The lower panel compares 
the binned values of (Aa/a) adc with the raw values from Fig. 6. 
Clearly, atmospheric dispersion effects (as modelled above) can not 
explain our results. Indeed, on average, the overall compression of 
the wavelength scale has an opposite effect on the low-z and high- 



10 We assumed the following observational parameters in our atmospheric 
dispersion model: temperature = 280 K, atmospheric pressure = 600 mbar 
and relative humidity = 10 per cent. 

11 In principle, one could find a reliable estimate of the seeing for each 
exposure by integrating the object profile in the spatial direction along the 
echelle orders. In general, we did not have access to the un-extracted, 2- 
dimensional object exposures and so could not perform such an analysis. 
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Figure 17. Correcting for potential atmospheric dispersion effects. The 
upper panel shows values of (Aa/a) ad derived from simulated spectra 
with our model for atmospheric dispersion applied. We have corrected the 
relevant raw values of Aa/a in the middle panel (solid squares) which 
also shows the unaffected raw values (open squares). The lower panel 
compares the binned values from the middle panel (solid circles) with 
those of the raw sample (open circles). The new overall weighted mean 
is Aa/a = (-0.46 ± 0.10) x 10" 5 . 



Table 6. Comparison of weighted mean values of Aa/a for the 
pre- and post-rotator systems. 
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Sample Pre-rotator 




Post-rotator 


7V abs Aa/a 1 


V ab5 


Aa/a 


Total sample (Fig. 20, middle panel) 






low-z 52 -0.60 ±0.16 


22 


-0.45 ± 0.20 


high-z 25 -0.39 ± 0.26 


29 


-0.83 ± 0.23 


Total 77 -0.54 ±0.14 


51 


-0.62 ± 0.15 


New sample only (Fig. 20, lower panel) 






low-z 22 -0.67 ± 0.25 


22 


-0.45 ± 0.20 


high-z 17 -0.77 ± 0.33 


17 


-0.50 ±0.31 


Total 39 -0.71 ± 0.20 


39 


-0.46 ± 0.17 



z systems. This is further borne out by the simulations in Section 
5.7.2 (see Fig. 23). 

In summary, it is clear that atmospheric dispersion effects can 
not explain the observed values of Aa/a. Assuming an effective 
seeing, our model of atmospheric dispersion provides a realistic 
upper limit to the potential effects on Aa/a. Additional track- 
ing/guiding errors will significantly reduce these effects by smear- 
ing the QSO light more evenly across the spectrograph slit. 



5.6.2 Evidence for atmospheric dispersion effects? 

If atmospheric dispersion significantly affects Aa/a then, from the 
discussion above, we expect to find a correlation between Aa/a 
and the zenith distance, £. We search for a such a correlation in 
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Figure 18. Is Aa/a correlated with zenith angle? (a) The estimates of the 
effect of atmospheric dispersion on Aa/a as a function of the mean zenith 
angle, f, for each observation. The Spearman rank correlation coefficient, 
r s , and associated probability, P(r s ), indicate a clear anti-correlation for 
the low-z systems but no clear correlation for the high-z systems, (b) The 
real raw values of Aa/a. No correlation is seen for either the low- or high- 
z systems, (c) A binning of the real values to clarify the lack of correlation 
with £. 
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Figure 20. Comparison of the pre- and post-rotator systems. The upper 
panel compares raw values of Aa/a for the pre-rotator (solid circles, 77 
points) and post-rotator (open circles, 51 points) samples. The weighted 
means are consistent: Aa/a = (—0.54 ± 0.14) X 10~ 5 and Aa/a = 
(—0.62 ± 0.15) X 10 — 5 respectively. The middle panel bins the results 
in the same fashion as for the main results in Fig. 6. We see no evidence 
for a significant systematic effect due to atmospheric dispersion. The lower 
panel compares the pre-rotator (solid squares, 39 points) and post-rotator 
(open squares, 39 points) subsamples of the new sample alone. Again, we 
see no evidence for the atmospheric dispersion effect. Table 6 presents the 
weighted mean values of Aa/a. 
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Figure 19. Is (Aa/ajjj correlated with Aa/a? The Spearman rank cor- 
relation coefficient and associated probability indicate no clear correlation 
for the low-z systems (solid squares). A weak correlation may exist in the 
high-z systems (open squares) but only at low significance. 



Fig. 18. In panel (a) we plot the estimates of (Aa/a) a d from the 
previous section versus the mean value of £ for each absorption sys- 
tem affected by atmospheric dispersion. The Spearman rank corre- 
lation coefficient, r B , and associated probability, P(r s ) 12 , indicate 
a clear anti-correlation for the low-z systems (solid squares). How- 
ever, we see no evidence for such a correlation in the real low-z 
values which are shown in panel (b). Panel (c) shows a binning of 
the results for clarity. This indicates that atmospheric dispersion has 
had no significant effect on the raw values of Aa/a. For the high- 



P(js) is the probability that |r s | could have been exceeded by chance. 



z systems, we also see no correlation between Aa/a and £ but, 
since the upper panel shows no clear correlation for the estimates 
of (Aa/a) a d, one cannot take this as evidence against atmospheric 
dispersion. One also expects a correlation between (Aa/a) a d and 
Aa/a if atmospheric dispersion effects are significant. However, 
in Fig. 19 we see no clear evidence for such a correlation for either 
the low- or high-z samples (solid and open squares respectively). 

We may also test for atmospheric dispersion effects by com- 
paring the value of Aa/a for the (77) pre- and (51) post-rotator 
values of Aa/a. In the top panel of Fig. 20 we identify these sub- 
samples and plot binned values in the middle panel. There is no 
discrepancy between the pre- and post-rotator samples. This is con- 
firmed in Table 6 where we compare the weighted mean values of 
Aa/a. Even for the low-z systems, which should be most sus- 
ceptible to atmospheric dispersion effects, we see no evidence for 
a discrepancy. We also consider the new sample separately in the 
lower panel of Fig. 20 where exactly half of the systems are pre- 
rotator observations. Again, in Table 6, the weighted mean values 
of Aa/a are consistent. 

To summarize this section, we emphasize that the QSO data 
themselves suggest that atmospheric dispersion had an insignificant 
effect on Aa/a. No correlation between Aa/a and the zenith an- 
gle of the QSO observations is observed (Fig. 18) and the pre- and 
post-rotator samples give completely consistent results (Fig. 20 and 
Table 6). 
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5.7 Line-removal tests 

A non-zero Aa/a manifests itself as a distinct pattern of line shifts 
(see Fig. 1). For a given absorption system, with a given set of fitted 
transitions, we can remove one or more transitions and still obtain 
a well-constrained value of Aa/a. In the absence of systematic er- 
rors associated with the transition removed or with the wavelength 
scale of the QSO spectra, Aa/a should not change systematically 
over all absorption systems. Removing a single transition or dis- 
tinct groups of transitions from our analysis therefore allows us 
to search for systematic errors without specific knowledge of their 
origin. Note that removing one or more transitions will result in a 
slightly modified estimate of the velocity structure, revising the es- 
timate of Aa/a for each system. However, this will have a random 
effect on Aa/a; the real question is, how robust are the values of 
Aa/a, averaged over the entire sample, to this line removal pro- 
cess? 

Below we construct three different line removal tests, each 
sensitive to different types of systematic errors. 



J. 7. 1 Single transition and species removal 

Removing a single transition (e.g. Mg II A2796) from the fits to the 
QSO data provides a direct means to search for systematic errors 
due to systematic line blending (Section 5.2.2), isotopic and hyper- 
fine structure effects (Section 5.5) and large errors in the measured 
laboratory wavelengths in Table 2. Removing entire species (e.g. all 
transitions of Si II) is also a test for isotopic and hyperfine structure 
effects. 

Given an absorption system, with a certain set of fitted transi- 
tions, we only remove transitions that can be removed. To clarify 
this, consider removing the Mg II A2796 transition. For the low-z 
Mg/Fe II systems, a well-constrained value of Aa / a can only be 
obtained when at least one Mg line is fitted since there is only a 
small difference between the q coefficients for the different (low- 
z) Fell transitions (see Fig. 1). Therefore, if the only anchor line 
present is Mg II A2796 then we can not remove it from the sys- 
tem. Similarly, we only remove the entire Mg II species if the Mg I 
A2852 line is fitted. 

We presented the results of removing single transitions and en- 
tire species for the previous samples in MOlb. In Fig. 21 we present 
line removal results for the total sample. The transition/species re- 
moved is given on the vertical axis together with the number of 
systems, n, where removal was possible. The left panel compares 
the weighted mean value of Aa/a before removing the transi- 
tion/species (dotted error bar) with that obtained afterwards (solid 
error bar) for these n systems. The right panel shows the effect of 
removing the transition/species on the weighted mean Aa/a for 
the sample as a whole (i.e. including those systems where the tran- 
sition in question could not be removed). Raw values of Aa/a are 
used in all cases (i.e. the 1 a errors returned from VPFIT are not 
increased to take account of extra scatter in high-z sample). 

Fig. 21 provides no evidence to suggest that systematic er- 
rors associated with any one transition or species have significantly 
affected Aa/a. This strongly suggests that isotopic/hyperfine sat- 
uration or evolution effects are not important for the transitions of 
interest. This is particularly important information for the majority 
of transitions for which the isotopic structures are not known (see 
Section 5.5). The line removal results also confirm that random and 
systematic blending with unidentified transitions has a negligible 
effect on Aa/a. 
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Figure 21. Single transition and species removal. The left panel compares 
the raw values of Aa/a before (dotted error bar) and after (solid error bar) 
line removal. The transitions or species removed are listed on the left to- 
gether with the number of systems, n, for which line removal was possible. 
The right panel shows the impact of line removal on the value of Aa /a for 
the entire sample. The vertical lines indicate the raw weighted mean value 
of Aa/a (dashed line) and 1 a error range (dotted lines). Note that some 
confusion may arise due to the occasional blending of the CrII and Znll 
A2062 lines: '(CrII + Znll) A2062' refers to cases where both transitions 
had to be removed simultaneously; 'CrII + Znll A2062' refers to cases 
when all CrII transitions were removed along with the blended Znll A2062 
line; a similar definition applies to 'Znll + CrII A2062'; 'Cm' and 'Znll' 
refer only to removal of the CrII and Znll species in cases where the CrII 
and Zn II A2062 lines were not blended. 



5.7.2 High-z compression test 

As previously noted, the arrangement of the q coefficients for the 
low-z Mg/Fe II systems (see Fig. 1) implies that a compression of 
the wavelength scale will systematically lead to Aa/a < 0. The 
dependence of Aa/a on the q coefficients is considerably more 
complicated for the high-z systems: a compression of the wave- 
length scale will have a different effect on each value of Aa/a 
depending on which transitions are fitted. That is, compression is 
not degenerate with Aa/a at high-z. If compression is responsible 
for the observed Aa/a < at low z, can evidence for it be found 
in the high-z absorption systems? 

We search for compression in the high-z sample by fitting 
combinations of transitions for which the q coefficients are ar- 
ranged in a similar way to those for the low-z Mg/Fe II systems. 
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Figure 22. The high-z compression test. A compression of the spectrum 
will systematically produce Act /a < in a high-z absorption system 
if only certain combinations of transitions (solid diagonal lines) are fitted 
simultaneously. This mimics the situation for the low-z Mg/Fe II absorp- 
tion systems and allows us to search for a compression-type systematic 
effect in the high-z systems. Combinations of transitions which will give 
Aa/a < for an expansion of the spectrum can be selected in a similar 
way (dashed diagonal lines). 



For example, consider a system where the following transitions are 
present in the QSO spectrum: Si II AA1526 & 1808, Fe II A A 1608 & 
1611 and CrII AA2056 & 2066. The arrangement of q coefficients 
for such a system is shown in Fig. 22. Fitting the Sill A1526 and 
Cr II transitions will yield Aq /a < if a compression of the spec- 
trum exists. Note that several different combinations could be used 
to the same effect, as indicated by solid diagonal lines. We treat 
such cases separately in our analysis since the lower wavelength 
transitions (i.e. Sill AA1526 & 1808 and Fell A1608) are of differ- 
ent q-type. The values of Aq / a and the 1 a errors in these separate 
cases will not be independent. Note also that we can form combi- 
nations of transitions which mimic Aa/a < for an expansion of 
the spectrum. 

From the high-z sample, 47 systems contribute 76 different, 
but not independent, combinations to the compression sample. Sim- 
ilarly, 59 systems contribute 126 combinations to the expansion 
sample. In Table 7 we calculate the raw values of the weighted and 
unweighted means for each sample both before and after removing 
transitions to form the compression and expansion combinations. 
Note that the number of compression combinations formed in a 
given absorption system is the number of times that system con- 
tributes to the pre-removal values. Table 7 shows no clear evidence 
for compression (or expansion) in the high-z sample. However, the 
weighted and unweighted means are not adequate statistics because 
each absorption system can contain many combinations. Further- 
more, the pre- and post-removal values and their 1 a errors are not 
independent of each other. Also of note are the large errors in the 
unweighted mean post-removal values, indicating a large scatter in 
the results and, therefore, that the weighted means are unreliable. 

To properly interpret results from the compression test, a com- 
parison must be made with detailed simulations of the QSO data. 
We synthesized each absorption system from the Voigt profile fit to 
the real QSO data, adding Gaussian noise appropriate to the mea- 
sured S/N. From each synthetic spectrum we constructed 3 different 
simulations: 
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Figure 23. Aa/a from simulations of the QSO spectra. The weighted 
mean from 20 simulations (solid circles) is compared with the real raw val- 
ues (dotted circles, shifted for clarity) in each panel. The rms from the simu- 
lations is represented by the error bar. The upper panel has no input value of 
Aa/a whereas Aa/a = —0.6 X 10~ 5 was introduced into the spectra for 
the middle panel. Note that we recover the input values and that the rms er- 
ror bars match the 1 cr errors in the real data. A compression was introduced 
into the spectra for the lower panel (equation 14) to allow interpretation of 
the high-z compression and (jr-type removal test results. 



(i) Aa/a = 0. No line shifts were introduced into the synthetic 
spectra. 

(ii) Aa/a = —0.6 x 10~ 5 . Each transition in the synthetic 
spectrum was shifted according to this value of Aa/a which cor- 
responds to the raw weighted mean value measured in the real QSO 
data. 

(iii) Compression. The synthetic spectrum was compressed ac- 
cording to 



A' = A - C*(A C - A) . 



(14) 



where A is the initial wavelength of a given pixel and A' is the 
new value obtained by compressing the spectrum about a central 
wavelength, A c , by a factor C. This mimics the compression effect 
of atmospheric dispersion, as described by equation 13. Thus, we 
chose A c = 5500 A as representative. In order to produce Aq /a « 
—0.6 x 10~ 5 in the low-z simulated spectra, we used C — 4 x 
10~ 6 . 

We ran 20 simulations of each absorption system, obtaining values 
of Aa/a for each of the 3 cases above. The results are compared 
with the real, raw sample in Fig. 23. Note that (i) we recover the in- 
put values of Aa/a for the Aa/a = and Aa/a = —0.6 x 10~ 5 
cases and (ii) the compression simulation shows a similar behaviour 
to that of the atmospheric dispersion calculation in Fig. 17, con- 
firming that simple distortions of the wavelength scale should have 
had an opposite overall effect on the high-z systems compared to 
the low-z systems. 

We have applied the high-z compression test to the simula- 
tions for each of the above 3 cases. The results are shown in the 
lower three panels of Fig. 24, binned for clarity. The error bars on 
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Table 7. Coarse results from the high-z compression test. Columns 4 and 5 show the raw values of the weighted 
and unweighted mean Aa/a (in units of 10~ 5 ) for the relevant absorption systems before removing the transitions. 
Each absorption system may contribute many times to these values, depending on how many compression/expansion 
combinations can be formed from the fitted transitions. Columns 6 and 7 give the weighted and unweighted mean after 
transitions are removed to produce the compression/expansion combinations. Note the large errors in the post-removal 
unweighted means. Other caveats of interpreting this table are discussed in the text. 
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Figure 24. Detailed results of the high-z compression test for the QSO data (upper panels) and simulations (lower 3 panels). The error bars for the simulations 
represent the rms from 20 synthetic spectra, (a) Compression sample. Each panel shows binned values of Aa/a before and after transitions are removed 
to form the 76 compression combinations. The bottom panel shows that the compression sample is sensitive to the synthetic compression introduced to the 
simulated spectra. No clear trend is seen in the real QSO data due to the large scatter in the post-removal values (see text for discussion), (b) Expansion sample. 
126 combinations contribute to the binned values of Aa/a shown. No significant difference is observed between the pre- and post-removal values in any of 
the simulations. The lower panel implies that the expansion sample is surprisingly insensitive to simple systematic errors. Significant scatter is also observed 
in the post-removal QSO results. 



the simulated pre- and post-removal values represent the rms from 
the 20 simulations of each absorption system. For the Aa /a = 
and Aa/a = —0.6 x 10~ 5 simulations, note that Aa/a does not 
change systematically after removing transitions to form the com- 
pression and expansion samples. As expected, the compression of 
the synthetic spectra causes the post-removal values of the com- 
pression sample (bottom panel of Fig. 24a) to shift towards more 
negative values. However, the compression simulation for the ex- 
pansion sample (bottom panel of Fig. 24b) reveals this sample to be 
surprisingly insensitive to the artificial compression of the spectra. 

The top panels of Fig. 24 compare the pre- and post-removal 
values of Aa/a for the real QSO data. The results are binned in 
the same way as the simulations but the error bars represent the 1 a 
error in the weighted mean of values within each bin. Of particu- 
lar note is the large scatter in the post-removal values compared to 
the size of the rms error bars in the simulations. The scatter is also 
large compared to the post-removal shift in Aa/a seen in the com- 
pression simulations. We explained the origin of the extra scatter at 



high-z in Section 4.2 and formed the fiducial sample by adding an 
additional random error to the systems likely to be most strongly 
affected. We expect this scatter to significantly increase as transi- 
tions are removed from the fit, as in the present case. The degree of 
extra scatter is, to some extent, represented by the large uncertain- 
ties in the unweighted means presented in Table 7. These values 
should therefore be far more reliable in such circumstances. 

It is therefore clear that the potential power of the compres- 
sion test (demonstrated in the bottom panel of Fig. 24a) is under- 
mined by the extra scatter in the high-z values. Unfortunately, one 
can draw no conclusions about potential systematic errors from the 
compression test without a significant increase in high-z sample 
size. 



5.7.3 High-z q-type removal 

If simple systematic errors are responsible for our measured non- 
zero Aa/a, it is surprising that Aa/a is so consistent between 
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Table 8. Coarse results of the high-z <?-type removal test. Columns 3 and 4 show the raw weighted and 
unweighted mean values of Aa/o (in units of 10~ 5 ) for the 26 relevant absorption systems before 
removing the transitions (the first three rows necessarily have the same value). Columns 5 and 6 give the 
weighted and unweighted mean after transitions of the specified g-type are removed. Note that the pre- 
and post-removal values of Act/a and the 1 a errors are not independent of each other. Also note the 
large errors in the post-removal unweighted means. A detailed comparison must refer to simulations of 
the QSO data (see Fig. 25). 
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Pre-removal 


Post-removal 






(Aa/a) w 


(Aa/a) 




(Aa/a) 


Anchors 


26 


-0.65 ±0.18 


-0.68 ±0.55 


-0.88 ±0.21 


-0.88 ± 0.46 


Positive-shifters 


26 


-0.65 ±0.18 


-0.68 ±0.55 


1.21 ±0.52 


0.25 ±1.39 


Negative-shifters 


26 


-0.65 ±0.18 


-0.68 ±0.55 


-0.63 ±0.34 


-1.89 ± 1.28 


Mediocre-shifters 


21 


-0.48 ±0.22 


-0.39 ±0.59 


-0.41 ±0.22 


-0.26 ± 0.72 



the low- and high-z samples. The lower panel of Fig. 23 explicitly 
demonstrates this. It would also be surprising if we were to take 
subsets of the QSO absorption lines, grouped according to the sign 
and magnitude of q, and find consistent values of Aa/a. We can 
apply such a test to the high-z sample only, since it contains transi- 
tions with several distinct types of q coefficient (see Fig. 1). 

We have applied this high-z g-type removal test to 26 absorp- 
tion systems. These contained at least one transition of each of 
the highly delineated g-types, i.e. anchors, positive- and negative- 
shifters. 21 of these systems also contain at least one mediocre- 
shifter. If a low-order distortion of the wavelength scale causes the 
observed non-zero Aa/a, removing all transitions of one g-type 
may result in a significant change in Aa/a. Table 8 compares the 
weighted and unweighted mean values of Aa/a before and after 
g-type removal. Although removing the anchors, mediocre-shifters 
and negative- shifters produces only small changes in Aa/a, re- 
moving the positive-shifters seems to cause to a significant change. 
However, two important caveats should be noted: (i) the pre- and 
post-removal values of Aa/a are not independent of each other. 
This also applies to the 1 a error bars, (ii) After removing the 
positive- and negative-shifters, the errors in the unweighted mean 
values are large, indicating a large scatter in Aa/a. Both these 
points make the results in Table 8 difficult to interpret. As in the 
previous section, the unweighted mean values should be far more 
reliable than the weighted means due to the extra scatter. This is 
confirmed below. 

As for the high-z compression test, comparison with simu- 
lations facilitates interpretation of the above results. We have ap- 
plied the g-type removal test to the 3 simulations described in 
Section 5.7.2. The results are shown in the lower three panels 
of Fig. 25, binned for clarity. The error bars on the simulated 
pre- and post-removal values represent the rms from 20 simula- 
tions of each absorption system. Firstly, note the consistency be- 
tween the pre- and post-removal values for the Aa/a = and 
Aa/a — — 0.6 x 10 simulations in panels (b), (c) and (d). Also, 
in the simulations of compressed spectra, we see that removing 
the anchor transitions [panel (b)] does not significantly affect the 
values of Aa/a. However, the effect is more significant when the 
positive- and negative-shifters are removed. Overall, we see Aa/a 
decrease when the positive-shifters are removed whereas removing 
the negative-shifters increases Aa/a. This illustrates the potential 
power of the g-type removal test in searching for low-order system- 
atic distortions of the wavelength scale. 

The top panels of Fig. 25b, c and d compare the pre- and post- 
removal results for the real QSO data. The results are binned in 
the same way as the simulations but the error bars represent the 



1 a error in the weighed mean of values within each bin. As for 
the high-z compression test results in Fig. 24, the results are con- 
fused by the significant scatter in the post-removal QSO values. 
We even observe extra scatter when removing the anchors from the 
real QSO data. Again, this extra scatter is expected. We also expect 
it to be worse here than for the compression test since (i) remov- 
ing the positive- or negative- shifters significantly reduces the sen- 
sitivity of the MM method to Aa/a, (ii) less absorption systems 
contain all three of the necessary g-types than systems contribut- 
ing to either the compression or expansion samples and (iii) many 
combinations of transitions can be used for each absorption system 
in the compression test. Also, by selecting systems which contain 
anchors, positive- and negative- shifters, we tend to select systems 
with transitions of very different line strengths (e.g. transitions of 
Alii, Cm and Znll; see Fig. 1). We explained how this leads to 
extra scatter in values of Aa/a in Section 4.2 and defined a 'high- 
contrast' sample of 22 systems at high-z. 18 of these systems are 
contained in the sample of 26 systems where g-type removal is pos- 
sible. The extra scatter is clearly apparent in Fig. 25a which shows 
the unbinned pre-removal values of Aa/a for these 26 systems. 

To summarize this section, simulations indicate that removing 
transitions, grouped according to the magnitude and sign of q, is a 
potentially powerful test for simple systematic errors in the high- 
z data. However, the extra scatter in the values of Aa/a for the 
relevant 26 absorption systems precludes any conclusion about the 
existence or magnitude of such systematic effects. 



6 DISCUSSION AND CONCLUSIONS 
6.1 Summary of new results 

We have used the MM method to analyse 3 samples of 
Keck/HIRES QSO spectra containing a total of 128 absorption sys- 
tems over the redshift range 0.2 < z a b s < 3.7. All 3 samples inde- 
pendently yield consistent, significantly non-zero values of Aa/a 
(Table 4, Fig. 6). Combining the samples gives a weighted mean 
Aa/a = (-0.574 ± 0.102) x 10~ 5 . For systems at z abs < 1.8, 
the predominant transitions arise from Mg II and Fe II whereas the 
higher z systems contain a variety of ions and transitions which 
have a diverse dependence on a. Therefore, if systematic effects 
were responsible for the observed non-zero Aa/a, one would ex- 
pect significantly different values at low- and high-z (Fig. 23). 
However, we find a similar Aa/a in both cases, suggesting that 
our results are robust. 

We identified a source of additional random scatter in the 
values of Aa/a for 22 (out of 54) high-z systems, i.e. those for 
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Figure 25. Detailed results of the high-z <jr-type removal test for the QSO data and simulations (lower 3 panels). The error bars for the simulations represent 
the rms from 20 synthetic spectra, (a) The upper panel shows the values of Aa/a for the 26 systems comprising the <?-type sample, i.e. those with at least 
one anchor, positive- and negative-shifter. These data are binned in the next panel. The lower three panels show the weighted mean values of Aa/a for the 
different simulations. These values and the binned QSO data are shown as dotted circles in (b), (c) and (d). (b) Anchors removed. The synthetic spectra with 
artificial compression (lower panel) indicate that removing the anchors should not affect Aa/a significantly, even when distortions of the wavelength scale 
are present. The QSO values also show little deviation after line removal but some extra scatter may exist, (c) and (d) Positive- and negative-shifters removed. 
Although the compression simulation reveals a significant change in Aa/a upon line-removal, the extra scatter in the post-removal QSO data precludes a firm 
conclusion about the presence of systematic effects (see text for discussion). 



which the fitted transitions have widely differing optical depths. In 
these systems, the velocity structure is primarily constrained by the 
low optical depth transitions. Any velocity components too weak 
to feature in low optical depth transitions may, nevertheless, sub- 
tly affect the fitted line positions of nearby velocity components in 
the high optical depth transitions. We have therefore increased the 
errors on Aa/a for these systems to match the observed scatter. 



Our most robust estimate for the overall weighted mean becomes 
Aa/a = (-0.543 ± 0.116) x 10~ 5 , representing 4.7 a statistical 
evidence for a smaller a in the QSO absorption systems (Fig. 8). 

If one assumes that Aa/a = at z a b s = then a constant 
increase in a with time is preferred by the data over a constant off- 
set from the laboratory value: a/a = (6.40 ±1.35) x 10 _16 yr _1 
(Fig. 9). However, a bootstrap analysis demonstrates that this pref- 
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erence has low significance. We find no evidence for dipolar vari- 
ations from the angular distribution of Aa/a (Fig. 10) and no 
evidence for spatial correlations in a over 0.2-13 Gpc comoving 
scales from the two point correlation function (Fig. 12). 

We have searched for possible instrumental and astrophysical 
systematic effects which could mimic the above evidence for vary- 
ing a. As we found in MOlb, the two most important of these are 
due to possible atmospheric dispersion effects and isotopic abun- 
dance variations. Sixty per cent of the absorption systems may be 
affected by the former. However, comparison of the affected and 
unaffected systems yields no evidence for these effects (Fig. 20). 
Modelling of atmospheric dispersion and its effect on the QSO 
spectra indicates that Aa/a in the low-z systems should be cor- 
related with the zenith angle of the QSO observations. We ob- 
serve no such correlation (Fig. 18). Atmospheric dispersion effects 
clearly can not explain our results. The effect of possible varia- 
tions in the isotopic ratios is more difficult to estimate. The iso- 
topic abundance trends in low metallicity stellar environments sug- 
gest that the 25 ' 26 Mg/ 24 Mg and 29 ' 30 Si/ 2S Si ratios will be signif- 
icantly lower in the QSO absorption clouds. This systematic ef- 
fect, if present, should have pushed Aa/a to positive values. Thus, 
when it is incorporated into the analysis, Aa/a becomes more neg- 
ative (Fig. 16). 

6.2 Comparison with other QSO absorption lines constraints 

The only other QSO absorption line method which directly con- 
strains Aa/a is the AD method outlined in Section 1.1.1. The MM 
method constraints in the present work are more than an order of 
magnitude more precise than the strongest AD method constraints 
(obtained in Murphy et al. 2001c). 

Comparison of H I 21 -cm and millimetre-band molecular ro- 
tational absorption frequencies constrains variations in y = a g p , 
where g p is the proton (/-factor (Drinkwater et al. 1998). This tech- 
nique suffers from an important systematic error which is diffi- 
cult to quantify: the millimetre-band and radio continuum emission 
generally originate from different regions of the background QSO, 
leading to possible line-of-sight velocity differences between the 
millimetre-band and Hi 21-cm absorption lines. Drinkwater et al. 
(1998) compared Galactic Hi 21-cm and millimetre-band absorp- 
tion lines to estimate the likely line-of-sight velocity difference, 
| Aw | = 1.2 km s -1 . This corresponds to a systematic effect for 
a single measurement of \Ay/y\ = 0.4 x 10~ 5 . This formed 
the dominant (1 a) uncertainty in our recent measurements of 
Ay/y for two absorption systems (Murphy et al. 200 Id): Ay/y = 
(-0.20±0.44) x 10~ 5 for2 abs = 0.2467 towards PKS 1413+135 
and Ay/y = (-0.16 ± 0.54) x 10~ 5 for z abs = 0.6847 towards 
TXS 0218+357. Carilli et al. (2000) argue that |Av| could be as 
large as 10-100 km s -1 on the basis of typical subkiloparsec mo- 
tions of the ISM. They derive the upper limit \ Ay/y\ < 1.7 x 10~ 5 
for the same two absorption systems. 

Comparison of optical and Hi 21-cm absorption frequen- 
cies constrains variations in x = a 2 g p m e /m p where m e /m p 
is the electron-to-proton mass ratio (Wolfe et al. 1976). Cowie & 
Songaila (1995) obtained Ax/x = (0.70 ± 0.55) x 10~ 5 from the 
Zabs = 1.7764 system towards Q1331+170. However, the effect of 
line-of-sight velocity differences has not been estimated for the op- 
tical/Hl comparison but is likely to be at least |Av| ~ 1.2 kms -1 , 
as found for the millimetre/H I technique. A statistical sample of 
Ax/x and Ay/y measurements is required to provide reliable con- 
straints and to properly quantify the typical line-of-sight velocity 
differences. This is currently limited by the small number of sys- 



tems known to exhibit both optical/millimetre-band and Hi 21-cm 
absorption (e.g. Curran et al. 2002). 

The above constraints on Ay/y and Ax/x can only be com- 
pared with MM constraints on Aa/a if one assumes the constancy 
of g p and g p m e /m p respectively. If one suspects variation in a 
then it seems wholly unjustified to assume the constancy of other 
fundamental parameters like g p and m e /m p . Some model depen- 
dent links between Aa/a, Ax / x, and Ay /y and have been studied 
recently within the paradigm of grand unification (e.g. Langacker, 
Segre & Strassler 2002; Dent & Fairbairn 2003; Dmitriev & Flam- 
baum 2003). Langacker et al. (2002) emphasize that the observa- 
tions should be used to constrain the theory and not vice versa. 

6.3 Comparison with non-QSO absorption lines constraints 

Several non-QSO absorption line constraints on Aa/a exist, 
falling into two distinct classes, 'local' and 'early universe'. We di- 
rect the reader to the recent review of Uzan (2003) for a summary 
of these methods. 

The local constraints include laboratory comparisons of 
atomic clocks (Prestage et al. 1995) and fountains (Sortais et al. 
2001), analysis of 149 Sm isotopic abundances from the Oklo nat- 
ural fission reactor in Gabon, Africa (Shlyakhter 1976; Damour 
& Dyson 1996), and meteoritic constraints on the long-lived (3- 
decay rate of 187 Re (Peebles & Dicke 1962; Dyson 1967, 1972). 
The strongest current constraints are, respectively, |Aa/a| = 
(0.2 ± 7.5) x 10~ 15 (r = 4.7yr, Marion et al. 2003), Aa/a = 
(-0.04 ± 0.15) x 10~ 7 (r = 1.8 Gyr, Fujii et al. 2000) and 
\Aa/a\ < 3 x 10~ 7 (r = 4.56 Gyr, Olive et al. 2002), where 
t is the relevant look-back time. However, none of these methods 
constrain a directly: the laboratory measurements assume constant 
nuclear magnetic moments (Marion et al. 2003), the Oklo bound 
can be weakened by allowing other interaction strengths and mass 
ratios to vary (Flambaum & Shuryak 2002; Olive et al. 2002) and 
Uzan (2003) notes that the meteoritic /3-decay limit assumes a con- 
stant weak interaction strength, a w . Moreover, as discussed in Sec- 
tion 4.3, one can not reliably compare these local limits with MM 
constraints on Aa/a without a detailed theory giving both tempo- 
ral and spatial variations of a. 

The early universe constraints come from analysis of the CMB 
anisotropies (Hannestad 1999; Kaplinghat, Scherrer & Turner 
1999) and light-element (e.g. D, He, Li) big bang nucleosynthe- 
sis (BBN) abundances (Kolb, Perry & Walker 1986). The strongest 
constraints from the CMB are at the Aa/a ~ 10~ 2 level if one 
considers the uncertainties in, and degeneracies with, the usual cos- 
mological parameters (i.e. f2 m , Q.\ etc.; e.g. Martins et al. 2003). 
However, Kujat & Scherrer (2000) and Battye, Crittenden & Waller 
(2001) note a crucial degeneracy between a and m e . Estimates 
based on the BBN abundance of 4 He suffer from a large uncertainty 
as to the electromagnetic contribution to the proton-neutron mass 
difference. The least model-dependent limits are those of Nollett 
& Lopez (2002), Aa/a = (3 ± 7) x 10~ 2 , who, like Bergstrom, 
Iguri & Rubinstein (1999), considered all light elements up to 7 Li, 
thereby avoiding this problem. However, Uzan (2003) notes that 
variations in a are degenerate with variations in a w , a s and G. 

6.4 Future MM method checks 

All the alternative methods discussed in the previous two sec- 
tions constrain a in (sometimes model-dependent) combination 
with other constants. Therefore, the only current avenue for con- 
fidently ruling out the present evidence for varying a is to obtain 
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independent MM constraints from QSO spectra. The present work 
has shown that no known systematic errors can mimic the effect of 
varying a in the Keck/HIRES spectra. However, if subtle, unknown 
instrumental effects cause the line-shifts we observe, the high qual- 
ity spectra now available from the VLT/UVES and Subaru/HDS 
will bear this out. If the present results are confirmed, iodine cell 
calibration techniques, similar to those used to identify extra-solar 
planets (e.g. Marcy & Butler 1992), could be applied to selected 
absorption systems to confidently rule out instrumental effects. 
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